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and  the  presence  of  a less  heat-resistant  fraction  was  developed 
and  tested.  The  model  was  developed  by  applying  systems 
analysis  to  population  dynamics  of  bacterial  spores.  A three- 
neck  flask  reactor  was  used  to  perform  isothermal  laboratory 
experiments  with  Bacillus  subtilis  spore  suspensions.  Results  of 
these  experiments  were  used  to  determine  various  constants  in 
the  mathematical  model.  The  model  was  then  tested  under 
nonisothermal  conditions,  and  showed  close  agreement  between 
predicted  and  measured  number  of  survivors  from  a lethal 
transient  temperature  process. 

Improved  biological  indicator  units  were  designed  and 
fabricated,  and  their  performance  compared  to  that  of 
traditional  units  used  in  commercial  practice.  The  fabricated 
units  offered  improved  heat  transfer  and  mechanical  properties 
that  allowed  the  spore  suspension  to  experience  the  actual 
thermal  history  of  the  product  with  added  protection  from 
pressure  changes. 

These  results  introduce  new  improved  methods  for  tne  design, 
calculation,  and  validation  of  thermal  sterilization  processes. 
Specifically,  they  establish  a procedure  by  which  the 
sterilization  effectiveness  of  a thermal  process  can  be  tested 
microbiological ly. 


INTRODUCTION 


Heat  inactivation  of  bacterial  spores  is  of  critical  impor- 
tance in  assuring  tHe  stability  and  safety  of  canned  foods  and 
parenteral  drugs.  Exposure  to  the  high  temperatures  of  commer- 
cial sterilization  processes  is  detrimental  to  the  nutritional 
and  organoleptic  characteristics  of  foods  and  to  the  stability 
of  drugs.  Consequently,  a heat  treatment  must  be  severe  enough 
to  accomplish  sufficient  inactivation  of  bacterial  spores,  while 
avoiding  excessive  exposure  of  the  product  to  high  temperature. 
Improved  methods  for  design  and  verification  of  thermal  proc- 
esses are  needed  to  ameliorate  the  precision  and  safety  of  these 
heat  treatments.  Development  of  Improved  methods  requires  an 
understanding  of  tne  engineering  sciences  of  heat  transfer  and 
the  reaction  kinetics  of  thermal  inactivation  of  bacterial 

The  design  of  thermal  sterilization  processes  has  been 
traditionally  based  on  the  assumption  that  thermal  inactivation 
of  bacterial  spores  can  be  modeled  as  a single  first  order 
reaction.  As  such,  it  may  be  described  by  a straight  line  when 
the  logarithm  of  the  number  of  surviving  spores  is  plotted 
against  time  of  exposure  to  a constant  lethal  temperature 
(survivor  curve).  However,  actual  survivor  curves  plotted  from 


laboratory  data  frequently  reveal  considerable  deviation  from  a 
single  straight  line  during  early  periods  of  exposure.  This  is 
a consequence  of  additional,  concomitant  reactions  such  as  heat 
activation,  injury,  and  preliminary  inactivation  (early  inacti- 
vation of  less  heat-resistant  fractions). 

Deviations  from  a straight  line  model  often  make  it  difficult 
to  test  or  validate  a thermal  process  microbiological ly.  Plate 
counts  from  a spore  suspension  subjected  to  a variable 
temperature  history  in  a commercial  process  often  do  not  agree 
Kith  predicted  results  from  the  single,  first-order  reaction 
model.  Also,  the  increasing  use  of  high  temperature-short  time 
sterilization  processes  with  aseptic  packaging  systems  requires 
better  kinetic  models  because  the  concomitant  reactions  may 
render  uncertain  the  degree  of  inactivation  achieved  during 
relatively  brief  exposures  to  lethal  temperatures.  Therefore, 
an  improved  kinetic  model  should  include  all  simultaneous, 
independent  reactions  in  order  to  predict  accurately  the  shape 
of  survivor  curves  obtained  from  laboratory  data. 

Bioindicator  units  (Bills)  offer  consioerable  potential  for 
microbiological  validation  of  thermal  processes  to  support  tne 
assurance  of  safe,  thermally-processed  products  by  providing 
direct  evidence  of  the  lethality  achieved.  Two  important 
requirements  for  adequate  performance  of  Bills  are  a)  the  spore 
suspension  held  within  a unit  must  experience  the  same 
temperature  history  as  the  coldest  place  in  the  product. 


and  b)  the  material  used  for  the  BID  should  be  Innocuous  to  the 
spores.  Therefore,  a bioindicator  unit  should  offer  the  least 
possible  resistance  to  heat  transfer  from  the  product  while 
remaining  well-insulated  from  heat  tnat  would  be  conducted  from 
the  can  wall-steam  interface,  and  it  should  be  nontoxic  to  the 
indicator  microorganism  contained  therein. 

The  need  for  more  complete  kinetic  models  and  improved 
bloindicator  units  for  the  biological  validation  of  tnermal 
sterilization  processes  led  to  the  objectives  for  this  work. 

The  objectives  were  a)  to  develop  and  test  an  improved  mathe- 
matical model  for  predicting  the  number  of  survivor  spores  in 
response  to  lethal  heat  treatment,  taking  into  account  the  con- 
current reactions  of  activation  and  injury  in  addition  to 
inactivation,  and  b)  to  design  and  fabricate  bioindicator  units 
with  improved  heat  transfer  characteristics  and  to  compare  their 
performance  with  that  of  traditional  units  used  in  commercial 
practice. 

Chronologically,  this  research  started  by  verifying  reports 
on  discrepancies  found  between  the  final  number  of  survivor 
spores  in  plastic  Bills  predicted  by  Bigelow's  model  and  measured 
by  microbiological  enumeration  for  a given  temperature  history. 
In  a second  stage,  the  relative  importance  of  heat  transfer  on 
the  performance  of  Bills  was  investigated  experimentally.  This 
involved  the  design,  fabrication  and  testing  of  Bills  witn 
different  heat  transfer  characteristics.  The  results  led  to  the 


conclusion  that  heat  transfer  cannot  completely  explain  the 
discrepancies  between  calculated  and  measured  numbers  of 
survivor  spores  after  a heat  sterilization  process  and  that  the 
conventional  model  utilized  to  describe  the  kinetics  of  bac- 
terial spore  inactivation  (8igelow's  kinetic  model)  is  not 
precise  enough  for  the  validation  of  thermal  sterilization 
processes. 

To  develop  systematically  an  improved  kinetic  model  of  spore 
inactivation,  consideration  was  given  to  the  fact  that  in 
addition  to  inactivation,  which  is  a transformation  of  viable 
spores  rendering  them  unable  to  grow  in  the  enumeration  medium, 
other  transformations  may  occur  which  are  not  considered  In  the 
models  utilized  in  practice  for  design  and  validation  of  thermal 
sterilization  processes.  Tne  main  factors  affecting  the  number 
of  survivor  spores  in  response  to  a thermal  treatment  have  been 
identified  as  activation,  injury,  the  existence  of  less  heat- 
resistant  fractions,  and  alternative  germination  mechanisms. 
Activation  is  a transformation  of  viable,  dormant  spores 
enabling  them  to  germinate  and  grow  in  an  enumeration  medium. 
Injury  is  a transformation  of  spores  that  modifies  their 
requirements  for  germination  and  growth. 

The  hypotheses  tested  in  this  research  were  a)  the  descrip- 
tion of  inactivation,  activation  and  injury  as  simultaneous  and 
Independent  first  order  reactions,  and  b)  the  existence  of  less 
heat-resistant  fractions  of  a spore  population,  and  c)  the 


existence  of  an  alternative  germination  system  whose  activation 
by  heat  allows  B.  subtilis  spores  to  grow  when  the  incubation 
temperature  is  32  deg.  C.  These  hypotheses  led  to  a conceptual 
model  of  a population  of  bacterial  spores  treated  at  lethal 
temperatures.  The  conceptual  model  was  analyzed  to  derive  a 
mathematical  model  which,  in  turn,  was  used  to  predict  the 
dynamic  responses  of  spore  populations  to  lethal  temperature 
regimes.  In  the  case  of  constant  lethal  temperature,  the 
mathematical  model  was  solved  to  obtain  formulae  for  the 
responses  (population  dynamics),  which  consisted  of  exponential 
functions  of  time  involving  parameters  corresponding  to  initial 
conditions  and  rate  constants. 

To  determine  values  for  parameters  of  the  model,  a series  of 
isothermal  experiments  was  performed  where  B.  subtilis  spores 
were  exposed  to  lethal  temperatures.  The  survivors  at  succes- 
sive times  were  enumerated  microbiologically,  and  the  experimen- 
tal responses  were  fitted  to  the  response  formulae  of  the  model 
using  a nonlinear  regression  technique.  Temperature  dependence 
of  rate  constants  was  found  to  be  well  described  by  Arrhenius' 
equation.  To  validate  the  model,  a dynamic  experiment  was 
performed  where  a spore  suspension  was  subjected  to  a variable 
temperature  history  that  also  was  used  to  generate  corresponding 
responses  of  the  model  by  simulation.  Hodel  generated  respon- 
ses were  compared  with  the  results  of  microbiological  enumera- 
tion of  survivor  spores  in  samples  taken  during  the  experiment. 


The  results  show  the  capability  of  the  new  model  to  describe 
well  the  dynamics  of  a population  of  bacterial  spores  exposed 
constant  and  transient  lethal  temperatures. 


LITERATURE  REVIEW 


Study  of  the  biological  validation  of  thermal  sterilization 
processes  using  bioindicators  has  important  physical  and 
biological  components.  Heat  transfer  plays  an  important  role  in 
the  ability  of  a bioindicator  unit  to  follow  the  temperature 
history  of  a product,  and  the  thermobacteriological  behavior  of 
a spore  suspension  used  as  the  indicator  determines  the  number 
of  survivors  to  lethal  heat  treatment.  This  review  deals  with 
information  pertinent  to  system  analysis  and  experimental  design 
corresponding  to  important  aspects  of  biological  validation  of 
thermal  sterilization  processes. 

In  the  first  section  of  this  chapter,  use  of  bloindicators  In 
the  validation  of  thermal  process  is  discussed,  with  emphasis  on 
basic  principles,  performance,  and  matters  of  practical 
importance.  The  state  of  the  art  in  validation  of  sterilization 
process  is  such  that  in  spite  of  the  potential  importance  of 
BIUs,  disagreement  between  the  number  of  survivors  determined 
from  the  traditional  model  and  microbiological  enumeration, 
often  by  several  orders  of  magnitude,  is  an  Important  restric- 
tion on  their  use.  The  most  extensive  section  presents  biologi- 
cal aspects  of  spore  population  dynamics  under  lethal  treatment 
temperatures.  The  complexity  of  biological  materials  requires 
careful  examination  of  information  available  on  the  species. 


transformations  and  interactions  involved.  This  information  is 
the  foundation  for  systematic  development  of  conceptual  and 
mathematical  models  that  allow  analyses  for  better  qualitative 
and  quantitative  understanding  of  the  behavior  of  spore  popula- 
tions. The  last  section  of  this  review  is  a brief  examination 
of  the  most  important  concepts  of  system  analysis  applied  to 
biological  systems. 

Use  of  Bioindicators  in  the  Validation  of  Thermal  Processes 


The  use  of  microbiological  methods  in  testing  the  effective- 
ness of  commercial  sterilization  of  food  and  pharmaceutical 
products  has  been  difficult  and,  at  best,  only  marginally 
effective.  Historical  efforts  made  use  of  inoculated  pact 
studies,  in  which  product  samples  were  inoculated  with  large 
concentrations  of  heat  resistant  spores  and  incubated  after 
processing  in  order  to  relate  the  fraction  of  product  spoiled  to 
the  number  of  surviving  spores.  Hore  recent  methods  use  small 
vials  containing  test  spore  suspensions  that  are  inserted  into 
product  containers  during  a sterilization  process.  The  vials 
are  subsequently  recovered,  and  the  survivors  are  enumerated 
microbiological ly  to  indicate  the  extent  of  product  steriliza- 
tion. However,  in  practice,  the  number  of  survivors  counted 
after  exposure  to  the  sterilization  process  seldom  agrees  with 


the  number  predicted  by  Bigelow's  kinetic  model  of  the  indicator 
spore  population  subjected  to  the  temperature  history  of  the 

Process  validation  has  been  a topic  of  particular  interest  in 
the  pharmaceutical  industry  over  the  last  ten  to  fifteen  years. 
"Process  Validation  is  establishing  documented  evidence  that  a 
process  does  what  it  purports  to  do"  (Agalloco  1986).  The  same 
author  reported  participating  in  validation  studies  of  steam 
sterilization  which  reduced  cycle  time  by  two-thirds  while  still 
maintaining  an  overkill  cycle,  thus  suggesting  that  very  large 
safety  margins  exist  in  commercial  sterilization  processes,  with 
consequent  waste  of  time  and  energy  and  negative  effects  on 
product  quality. 

Biological  indicators  are  important  tools  in  tne  validation 
of  sterilization  processes.  These  are  standardized  preparations 
of  specific  microorganisms  of  known  resistance  to  a steriliza- 
tion process  agent.  In  the  purest  form  of  bioindication, 
microorganisms  are  added  directly  to  the  product  being  sterili- 
zed. When  this  approach  is  not  practical,  an  alternative  is  to 
add  the  organisms  to  carrier  materials,  or  containers,  known  as 
bioindicator  units  (BlUs).  In  appropriate  applications  of  this 
type  of  indicator.  Bills  should  have  a time  constant  of  heating 
and  cooling  that  allows  the  contents  to  follow  closely  the 
thermal  history  of  the  product  (Joslyn,  1983). 


Ill 

There  are  two  categories  of  tests  by  which  the  efficacy  of 
sterilization  can  be  assessed:  either  the  product,  or  a statis- 
tically significant  sample  of  it,  is  directly  tested  for  the 
presence  of  viable  microorganisms,  or  the  time-temperature 
conditions  of  the  sterilization  process  are  monitored  and 
analyzed  in  terms  of  existing  knowledge  about  the  behavior  of 
microorganisms  under  the  stress  conditions  of  sterilization.  In 
this  regard.  Bills  have  the  advantage  over  other  monitors  of 
responding  to  all  the  component  factors  of  a process.  The 
essential  requirement  for  effective  performance  of  a BIU  is  that 
its  behavior  under  stress  shall  be  predictable  (Quesnel,  1984). 
Pflug  and  Oudlaug  (1986)  published  a review  on  the  use  of  Bills 
in  the  pharmaceutical  and  medical  industry.  A biological 
indicator  system  is  defined  to  consist  of  microorganisms, 
usually  bacterial  spores,  procedures  and  equipment  for  preparing 
the  Bill,  and  the  measurement  procedure  for  determining  tne  level 
of  sterilization  achieved  at  a specific  location  in  the  product. 

8IUs  that  use  living  bacterial  spores  as  sensing  units  cannot 
be  considered  a calibration  standard.  From  a standpoint  of 
metrology.  Bills  must  be  calibrated  in  a laboratory  using 
systems  that  are  based  on  physical  standards.  Pflug  and  Oudlaug 
(1986,  pages  244-245)  also  mention  that  there  is  general 
understanding  and  agreement  among  scientists  working  on 
sterilization  in  the  pharmaceutical  industry  and  in  regulatory 
agencies  that  "the  design  of  a sterilization  process  for  a 


product  that  will  be  marketed  as  sterile  must  be  based  on  the 
number  of  microorganisms  on  or  in  the  product  (bioburden)."  the 
probability  of  a nonsterile  unit  s(PNSU)  on  the  order  of  10“6  is 
widely  accepted  as  a standard  target  for  most  terminally- 
sterilized  products.  The  conclusion  is  that  Bills,  with  their 
inherent  variability,  are  reasonably  accurate  and  reproducible 
in  measuring  a sterilization  effect  when  properly  used. 

Hyers  and  Chraf  (1980)  published  a review  on  the  basic 
principles  and  applications  of  bioindicators.  Bioindicator 
systems  should  be  prepared  to  offer  a realistic  challenge  to  the 
process.  One  approach  is  to  set  the  concentration  of  the 
indicator  organism  to  the  bioburden  of  the  material  to  be 
sterilized.  One  may  exceed  the  bioburden  concentration  to 
provide  a safety  factor.  Care  should  be  taken  not  to  create 
situations  that  would  add  resistance  to  penetration  of  heat  from 
the  sterilizing  medium.  Bruch  (1974)  published  an  analysis  of 
the  probabilities  of  survivors  which  showed  the  advantage  of 
BlUs  over  traditional  sampling  methods.  The  review  was  motiva- 
ted by  a series  of  septicemias  originated  from  contaminated, 
large  volume,  parenteral  solutions.  The  main  conclusion  of 
8ruch's  study  is  that  the  simplest  way  to  achieve  the  desired 
probability  of  10"6  (practically  impossible  by  increasing  sample 
size)  is  through  the  use  of  Bills. 

Plastic  Bills  are  described  and  discussed  by  several  authors 
(Jones  et  al.,  1980,  Pflug  et  al.,  1980  and  Pflug,  1982).  Jones 


1.(1980)  describe 


fill  weight  on  the  F-value  delivered  to  green  beans  using  a 
sterilmatic  retort.  Pflug  et  al.  (1980)  compared  differences  in 
the  values  of  lethalities  calculated  by  the  general  method  (Ball 
and  Olson,  1957)  and  determined  by  microbiological  enumeration 
of  survivor  spores  of  8.  stearothermophilus.  They  found  them  to 
be  far  larger  than  variations  due  to  experimental  error  ( 1.2 
min.  vs  0.03  to  0.06).  Average  differences  between  calculated 
and  measured  lethalities  were  sometimes  as  big  as  2.4  min.  The 
authors  said  the  conclusions  that  could  be  made  were  limited  due 
to  the  complexity  of  the  overall  measurement  problem. 

The  United  States  Pharmacopoeia  (USP)  has  supported  and 
encouraged  studies  aimed  at  setting  standards  for  BIUs.  Cooper 
(1985)  introduced  a series  of  reports  and  monographs  from  the 
USP  that  have  been  published.  Several  of  the  reports  reveal  the 
need  to  improve  the  quality  of  81Us  commercially  available 
(Hayernik,  1972;  Hacek,  1972).  Smith  et  al.  (1982)  published  a 
report  on  three  different  commercially  available  biological 
Indicator  spore  strips.  Only  one  type  met  the  manufacturer's 
specifications.  Gil  Its  (1975)  tested  six  commercially  available 
bioindicator  systems.  He  found  that  not  all  of  them  were 
properly  labelled  and  that  the  claims  of  some  manufacturers  were 
not  sustained.  Vawger's  (1978)  paper  on  bacteriological 
evaluation  of  thermal  process  design  presents  the  approach  of 
the  food-processing  industry,  which  is  based  on  the  inoculated 


pack  system  and  on  the  count  reduction  system.  The  product- 
fractional  exposure  method  for  validation  of  sterilization 
processes  was  studied  by  Caputo  et  al.,  (1979).  In  this  method, 
nonsterile  finished  product  is  exposed  in  a production  vessel  to 
various  fractions  of  a proposed  or  existing  sterilization  cycle 
time,  and  the  surviving  bioburden  is  measured  for  each  fraction 
of  the  cycle  to  determine  the  time  necessary  to  achieve  sterili- 
zation. The  use  of  B.  stearothermophilus  as  a biological 
indicator  was  investigated  by  Henning  et  al.  (1973).  The  effect 
of  storage  conditions  on  Bids  using  8,  stearothermophilus  was 
explored  by  Belch  et  al.  (1979)  over  a 2.5  year  period  that 
demonstrated  the  decline  in  resistance  over  time  of  several 
commercially  available  BIUs. 

Bacterial  Spore  Inactivation 

All  species  of  bacterial  sporeformers  are  in  the  family 
Bacillaceae.  divided  into  the  five  genera:  Bacillus.  Clostri- 
dium, Sporosarcina.  Oesulfotomaculum.  and  Sporolactobacillus. 
Since  the  simultaneous  discovery  of  bacterial  spores  by  Koch  and 
Cohn  in  1876,  a considerable  amount  of  research  on  spores  has 
been  done.  Studies  on  thermal  inactivation  of  bacterial  spores 
are  historically  related  to  studies  on  disinfection.  In  1897, 
Kronig  and  Paul  published  an  extensive  paper  on  disinfection 
that  established  some  of  the  basic  criteria  and  procedures  for 


experimental  study  of  disinfection:  the  effects  of  time, 
concentration  and  type  of  disinfectant,  temperature,  and  type  of 
test  bacteria  on  the  number  of  survivors  were  carefully  examin- 
ed, and  several  rules  for  testing  the  efficacy  of  disinfectants 
were  proposed.  The  authors  were  aware  of  the  work  of  Arrhenius, 
and  they  recommended  the  use  of  physical-chemistry  in  the  study 
of  disinfection. 

In  1910  Chick  published  a paper  on  disinfection,  where  hot 
water  was  used  as  a disinfectant.  The  author  reviewed  the 
results  published  by  Kronig  and  Paul  (1897)  and  showed  tnat 
disinfection  proceeded  in  an  orderly  manner  and  that  the  rate 
diminished  as  the  number  of  survivors  became  less.  Mention  was 
made  that,  in  previous  work  with  mercuric  chloride  and  phenol, 
the  concentration  of  spores  remaining  alive  varied  logarithmi- 
cally with  time  and  that  the  number  of  spores  destroyed  per  unit 
time  was  proportional  to  the  concentration  of  viable  spores  at 
that  moment.  Analysis  of  the  data  published  by  Kronig  and  Paul 
(1897)  showed  the  same  relationship.  Experiments  were  conducted 
to  determine  first-order  rate  constants,  and  an  attempt  was  made 
to  use  Arrhenius'  equation  between  45  and  55  deg.  C.  Some 
deviations  from  logarithmic  linearity  were  found.  Similitude  of 
the  denaturalization  of  proteins  under  heat  treatment  was  noted, 
together  with  very  high  values  of  the  activation  energy.  The 
conclusion  uas  that  disinfection,  whether  by  disinfectants  or  by 


considered  analogous 


chemical  reaction. 


rate  of  which  is  controlled  by  conditions  of  temperature  and 
concentrations  of  bacteria  and  disinfectant. 

Viljoen  (1924)  studied  the  protective  effect  of  sodium 
chloride  on  bacterial  spores  heated  in  pea  liquor.  This  work 
related  directly  to  the  canning  industry,  dealing  with  the 
effect  of  product  composition  on  spore  inactivation.  The 
protective  Influence  was  apparent  to  sodium  chloride  concentra- 
tions of  3 to  3.5X.  Higher  concentrations  had  no  protective 
effect  and  even  produced  a killing  effect.  The  logarithmic 
character  of  spore  inactivations  also  was  noted,  and  the  results 
were  presented  in  semi-logarithmic  form.  The  importance  of 
environment  composition  was  also  studied  by  Watkins  and  Winslow 
(1932),  and  survivor  curves  were  found  to  be  logarithmic  in 
nature.  The  effect  of  alkalinity  on  rate  constants  was  such 
that  doubling  the  concentration  of  alkali  increased  the  value  of 
rate  constants  about  tenfold,  and  increasing  the  temperature 
from  50  to  55  deg.  C increased  the  value  of  rate  constants 
twenty  to  thirtyfold.  The  effects  on  rate  constants  of  age  and 
initial  concentration  of  the  inoculum  were  studied,  and  ageing 
was  found  to  decrease  while  concentration  increased  the  values 
of  rate  constants. 

The  importance  of  predicting  inactivation  of  bacterial  spores 
for  the  canning  industry  led  to  the  development  of  empirical 
methods.  The  thermal  death  point  was 


defined  as  the 


interval  at  a given  temperature  necessary  to  destroy  a specified 
initial  concentration  of  spores  in  a medium  of  known  pH,  and  a 
method  was  proposed  for  determining  the  thermal  death  point 
(Bigelow  and  Esty,  1920).  Later,  the  logarithmic  nature  of 
thermal  death-time  curves  was  finally  taken  into  consideration, 
and  tne  use  of  semi-logarithmic  graphs  was  proposed  for 
determining  the  time  necessary  to  destroy  a certain  fraction  of 
the  initial  viable  spores  (Bigelow,  1921). 

The  General  Method  (Bigelow  et  al.,  1920)  for  design  and 
evaluation  of  thermal  processing  foods  is  a graphical  integra- 
tion procedure  for  determining  the  lethal  effects  of  different 
time-temperature  regimes  at  the  coolest  spot  in  a container. 
Computers  make  this  method  easy  to  use  today,  but  it  was  tedious 
and  time-consuming  when  first  introduced.  Ball  (quoted  in 
Stumbo,  1948)  developed  a mathematical  method  for  calculations 
pertaining  to  thermal  processes  that  is  frequently  used  today 
(Formula  Method).  Ball  was  a strong  proponent  of  using 
mathematical  tools  in  practice;  his  formula  method  is  based  on  a 
hyperbolic  approximation  to  the  heat  penetration  curve  in  a 
finite  cylinder  (cylindrical  can).  For  the  same  purpose,  Stumbo 
(1948)  developed  a mathematical  method  that  considers  all  points 
in  the  container.  Tnis  was  an  improvement  over  Ball's  formula 
method,  which  considered  only  the  can  center  as  the  slowest- 
heating  point.  Therefore,  by  1950,  several  methods  were 
available  that  allowed  calculations  for  thermal  sterilization 


processes.  This  subject  was  thoroughly  analyzed  and  discussed 


in  two  books  that  may  be  considered  classics;  they  provide  ample 
information  on  the  theoretical  and  practical  aspects  of  thermal 
process  design  (Ball  and  Olson,  1957,  and  Stumbo,  1965).  Ball 
(1938)  also  discussed  the  relationship  between  thermal  processes 
and  the  organoleptic  characteristics  of  canned  products,  and  he 
recommended  high-temperature  short-time  treatments  of  natural 
and  forced  convection-heating  foods.  The  above  mentioned 
works  form  the  basis  for  modern  design  and  evaluation  of  thermal 
sterilization  processes,  and  tney  rely  on  pseudo  first-order 
kinetics  to  describe  inactivation  of  bacterial  spores. 


The  traditional  view  of  bacterial  spore  inactivation  at  a 
constant  lethal  temperature  may  be  represented  by  the  model  in 
Figure  1 composed  of  a store,  N,  of  viable  spores,  an 
inactivation  transformation,  0,  and  a sink  (triangle)  for  dead 
or  inactivated  spores.  In  Figure  1,  the  variables  denoted  by 
identify  flows  related  to  these  processes.  An  analysis  of  the 
model  is  now  presented. 

The  description  of  the  store  is 
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Conceptual  model  corresponding  to  the  single 
transformation  approach  to  spore  inactivation. 
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Inactivation  as  an  algebraic  transformation  is 

''O'KN  (2) 

where  K is  the  rate  constant  and  t is  time.  The  mode]  has  a 
single  node,  at  which  the  sum  of  incident  flows  equals  zero: 

-YD  - \ • 0;  or  1 N ■ -To  (3) 

Substitution  of  equations  1 and  2 in  3 yields  the  overall  model 
of  inactivation  as  analogous  to  a first  order  reaction: 


■(4) 


When  the  initial  number  of  viable  spores  is  No  the  solution 
of  (4)  is  as  follows: 


N = No  exp(-K  t)  (5) 

The  logarithmic  form  of  the  solution  often  is  used, 
ln(N/No)  = -It  t,  or 

log(No/N)  = t/D  (6) 

where  D,  decimal  reduction  time,  is  the  time  interval  required 
to  reduce  the  population  of  viable  spores  to  one-tenth  of  its 
former  value  ( note  that  0 ■ ln(10)/K).  Figure  2 is  a 


Figure  2.  Semi-logarithmic  survivor  curve.  Decimal 
reduction  time,  0,  is  the  time  required  by 


magnitude. 
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semilogarithmic  plot  of  the  inactivation  of  a bacterial  spore 
population,  the  thermal  death  time  (TOT)  curve,  where  0 is 
indicated  to  be  the  Interval  required  for  the  TOT  curve  to 
advance  one  log  cycle. 

In  practice,  thermal  processes  are  often  described  by  their 
lethalities.  For  a constant  treatment  temperature,  lethality  F 
is  defined  to  be  the  time  required  for  a specified  reduction  of 
the  population  of  viable  spores.  From  Equation  6,  lethality  at 
a constant  temperature  is  expressed  by: 

F = D log(No/N)  (7) 

Since  0 is  temperature  dependent,  F is  expressed,  by  convention, 
for  a reference  temperature  of  121.1  deg.  C (250  deg.  F): 

f121.1  ■ 0121.1  log(No/N)  (8) 

For  instance,  if  the  value  of  the  decimal  reduction  time  at 
121.1  deg.  C is  3 min  and  the  desired  reduction  in  the  number  of 
viable  spores  is  9 orders  of  magnitude,  the  lethality  of  the  is 

The  empirical  expression  used  in  practice  to  describe  the 
variation  of  0 with  temperature  is 


■(9) 


log(0n/0T2)  . (T2  - TiJ/a 
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where  T i and  Tj  are  two  different  lethal  temperatures  and  z is 
the  temperature  increment  required  to  reduce  the  value  of  the  0 
by  one  order  of  magnitude,  as  shown  in  Figure  3. 

For  a variable  temperature  process,  lethality  Is  determined 
by  dividing  the  thermal  history  into  small  time  intervals 
wherein  the  temperature  is  considered  constant.  The  effect  of 
temperature  during  an  interval  (t2  - tj)  on  the  number  of 
survivor  spores  follows  from  equation  6: 

log(N!/N2)  = (t2  - ti)/Dt2 (10) 

If  the  time  increment  is  of  differential  magnitude,  the  total 
lethal  effect  is  obtained  by  adding  the  contributions  of  the 
differential  treatments: 

log(No/li)  = , “ 

' (11) 

D(T) 

From  the  description  of  temperature  dependence  of  0 (equation  9) 
and  using  the  121.1  deg  C (250  deg.  F)  reference  temperature, 
the  value  of  1/D(T)  is 

(T2-121.1)/z 

1/D(T)  * (1/Otfi.!)  * 10  ..  (12) 


Figure  3.  Thermal  destruction 
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The  terms  defined  — D,  z and  F—  are  commonly  used  to  define 
the  heat  resistance  of  bacterial  spores  and  the  lethal  effect  of 
sterilization  processes  on  them. 

Figure  4 presents  an  example  of  the  complex  behavior  that  may 
be  observed  during  thermal  inactivation  of  bacterial  spores;  it 
includes  an  initial  decline  in  the  number  of  survivors  due  to 
preliminary  inactivation,  a subsequent  rise  in  the  number  of 
survivors,  and,  finally,  the  predominant  inactivation  phase. 

The  straight  line  (I)  corresponds  to  the  simple,  traditional  model 
(Equation  6)  where  inactivation  is  assumed  to  be  the  only 
significant  reaction.  It  also  corresponds  to  the  related 
practice  of  defining  the  kinetics  of  spore  inactivation  at 
constant  lethal  temperature  by  fitting  Equation  6 to  only  two 
points:  at  time  equal  to  zero  and  at  the  time  when  the  popu- 
lation has  been  reduced  to  one  viable  individual  per  unit  volume. 
In  this  case,  the  most  probable  number  (MPN)  of  survivors  tech- 
nique is  used  to  enumerate  survivor  spores.  Use  of  the  simple 
straight  line  model  (Equation  6)  may  lead  to  large  errors  in  the 
estimation  of  bacterial  spores  that  will  survive  a lethal  heat 
treatment  or  the  heat  treatment  neccessary  to  render  a stable 
and  safe  product. 

The  modeling  of  complex  systems  of  irreversible  first-order 
chemical  reactions  was  investigated  by  Lee  (1978).  A lumped 
structure  allowed  for  the  synthesis  of  kinetically  consistent 
structures  of  reaction.  That  paper  dealt  with  organic  chemical 


Figure  4.  Deviations  from  logarithmic  linearity 

during  inactivation  of  bacterial  spores. 

I corresponds  to  the  logarithmic  behavior,  II  to 
preliminary  inaccivation  of  less  heat-resistant 
fractions.  III  is  the  region  dominated  by 
activation  and  IV  corresponds  to  predominant 
inactivation. 
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reactions.  Prokop  and  Humphrey  (1970)  presented  a review  paper 
on  models  of  disinfection  kinetics.  The  proposed  models  share  a 
consecutive-reaction  mechanism.  Shull  et  al.,  (1963)  also 
proposed  a consecutive  reaction  mechanism  for  spore 
inactivation,  postulating  the  existence  of  activated,  interme- 
diate, and  inactivated  states.  As  a consequence,  activation  and 
inactivation  would  be  dependent,  i.e.,  the  spores  must  be 
activated  first  and  then  inactivated.  No  biochemical  basis  was 
presented  for  this  assumption. 

Hore  recently,  Kalinina  and  Motina  (1984)  developed  a kinetic 
model  for  B.  stearothermophilus.  that  also  assumed  sequential 
reactions.  Arabshahi  and  Lund  (1985)  published  a paper  on 
procedures  for  properly  calculating  kinetic  parameters  from 
experimental  data.  Toda  and  Aiba  (1966)  investigated  the 
thermal  characteristics  of  spore  clumps.  The  Increase  in  life- 
span of  a clump  was  proportional  to  the  logarithm  of  the  number 
of  single  spores  within  the  clump.  Aiba  and  Toda  (1965) 
analyzed  thermal  death  rates  of  bacterial  spores.  They  found 
that  a normal  distribution  of  individual  life-spans  produced  a 
thermal  death  rate  of  a population  of  spores  in  the  same  form  as 
Equation  6.  The  life-span  of  a spore  was  defined  as  the  length 
of  time  the  spore  remained  viable  in  a given  temperature.  Each 
spore  was  assumed  to  have  a specific  life-span,  and  a spore 
suspension  was  assumed  to  have  some  distribution  of  life-spans. 
The  authors  recommended  supplementing  the  first-order  reaction 
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rate  model  with  probability  theory.  Berry  and  Bradshaw  (1982) 
compared  the  lethalities  of  thermal  sterilization  processes 
obtained  from  thermal  histories  using  Bigelow's  model  and  from 
microbiological  enumeration  of  survivors.  They  found 
differences  as  large  as  7.2  min.  (the  calculated  value  was 
larger  than  the  microbiological  measurement),  and  they  concluded 
that  the  two  sterilization  values  are  not  directly  comparable. 

Equipment  for  the  experimental  study  of  inactivation  kinetics 

The  need  for  precise  control  of  exposure  to  lethal  tempera- 
tures in  the  measurement  of  thermal  inactivation  kinetics  has 
led  to  considerable  research  with  alternative  laboratory 
equipment  and  methodology.  Srimani  and  Loncin  (1980)  used  a 
three-neck  flask  reactor  to  measure  survivor  spores  for  a 
treatment  temperature  of  100  deg  C.  The  survival  curve  for  8^ 
stearothermophilus  showed  a maximum  after  4 hr.  of  exposure  (one 
log.  cycle  increment),  when  ca.  10*  of  the  spores  found  by 
microscopic  count  were  activated.  8urton  et  al.,  (1977)  studied 
the  thermal  death  kinetics  of  B.  stearothermophilus  at  ultrahign 
temperatures  (UHT).  The  capillary  tube  method  was  compared  with 
a direct  heating  UHT  plant.  Perkin  et  al.,  (1977)  analyzed  the 
heat  transfer  of  the  capillary  tube  method  to  determine  possible 
correction  factors  for  UHT  treatments.  Oavies  et  al.,  (1977) 
found  that  capillary  tubes  could  be  used  up  to  132.5  deg.  C 


Arrhenius’  kinetics.  A continuous  flow 


before  departing  from 
apparatus  for  the  determination  of  heat-inactivation  kinetics 
was  published  by  Oquendo  et  al.,  (1975)  and  Srimani  et  al., 
(1980).  Kooiman  (1974)  proposed  the  use  of  a screw-cap  tube 
technique  where  an  aliquot  of  spore  suspension  was  injected  into 
a heating  medium  at  the  treatment  temperature.  The  volume  of 
the  heating  medium  was  99  times  larger  than  that  of  the  spore 
suspension;  this  permitted  very  fast  temperature  rise  (the 
change  in  the  reactor  internal  temperature  is  of  the  order  of  IX 
and  process  temperature  is  regained  in  a few  seconds).  Nallidis 
and  Scholefield  (1985)  developed  a solid  heating  block  system 
for  determining  heat  resistance  of  spores.  Erdtsleck  and  Beumer 
(1976)  proposed  a TDT-pouch  method.  David  and  Shoemaker  (1983) 
developed  a computer-controlled  reactor  using  an  air-operated 
piston  to  move  into  exposure  a known  volume  of  sample. 

The  use  of  differential  scanning  calorimetry  in  studies  of 
bacterial  thermal  death  has  been  investigated  by  several  authors 
(Grieme  and  Barbano,  1983,  and  Hiles  et  al.,  1986).  The 
thermograms  of  highly  heat-resistant  spores  are  different  from 
those  of  less  heat-resistant  spores.  Differential  scanning 
calorimetry  may  become  an  important  research  tool,  since  the 
measured  changes  in  heat  capacity  are  related  to  changes  in  the 
degrees  of  freedom  of  molecules  that  may  be  involved  in 
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rearrangements  of  vital  spore  structures,  (tiles  et  al.,  (1986) 
developed  equations  for  differential  scanning  calorimetry  to 
predict  survivorsliip  when  the  heating  rate  is  a constant. 

Physioloov  of  Bacterial  Spores 

Hodern  design  and  evaluation  of  thermal  sterilization 
processes  is  based  on  the  simple  pseudo  first-order  kinetics  of 
the  inactivation  of  bacterial  spores.  In  order  to  achieve  a 
better  understanding  of  the  various  reactions  that  occur  during 
the  exposure  of  bacterial  spores  to  lethal  temperatures, 
relevant  aspects  of  bacterial  spore  physiology  will  be  reviewed 
in  more  detail  in  the  following  subsections. 

Murrell  (1969)  presented  a comprehensive  review  of  the 
structure  and  chemical  composition  of  bacterial  spores  that  is 
frequently  mentioned  in  the  literature.  Reference  was  made  to 
the  fact  that  bacterial  spores  are  complex  in  structure  and 
composition,  differing  from  vegetative  cells  by  the  presence  of 
various  additional  layers,  unique  substances  such  as  diplcolinic 
acid  (DPA) , and  quantitative  differences  such  as  a 1000  times 
higher  calcium  concentration.  The  preparation  of  spores  for 
analysis  was  discussed  in  detail.  Bacterial  spores  contain  about 
10Z  nitrogen,  5*  carbohydrate,  15*  calcium  dipicolinate,  and  8* 
ash.  Spores  of  all  species  have  the  same  basic  structure:  A 
central  core  or  protoplast  surrounded  by  the  plasma  membrane. 


plasma 


germ  cell  wall,  cortex,  coats  and  exosporium.  The  core, 
membrane  and  germ  cell  wall  constitute  a condensed  cell, 
contained  and  protected  by  the  outer  structures  (Warth,  1978a). 
The  exosporium  may  be  tightly  or  loosely  fitting.  In  cases 
where  information  was  available,  the  exosporium  was  found  to  be 
mostly  protein  with  significant  amounts  of  carbohydrate. 

Spore  coats  show  an  interesting  variety  in  both  appearance 
and  complexity.  Three  main  types  of  layers  usually  can  be 
distinguished.  The  most  distinctive  is  the  middle  layer,  which 
shows  a characteristic  laminar  pattern.  Some  species  have 
heavily  rigged  and  ornate  coats,  while  others  have  simple,  thin 
coats.  Coats  consist  largely  of  structural  protein.  Smaller 
amounts  of  carbohydrates  and  lipids  are  generally  found,  and  in 
some  species  large  amounts  of  phosphorus.  Coats  comprise  30  to 
60*  of  the  dry  weight  and  40  to  80*  of  the  proteins  of  the 
spore.  The  major  extractable  protein  appears  to  be  similar  in 
distinct  species. 

Construction  of  spore  coats  is  a major  activity  of  a mother 
cell  during  spore  formation.  Two  main  stages  may  be  distin- 
guished: first,  biosynthesis  of  polypeptides  and,  second,  the 
aggregation,  modification,  addition  and  re-arrangement  involved 
to  assemble  structures.  Spore  coats  first  appear  after  the 
forespore  has  contracted  and  concomitantly  with  cortex  formation 
and  the  beginning  of  refractility.  Completion  of  the  coats  is  a 
late  event,  occurring  at  the  same  time  as  full  refractility  and 


heat  resistance.  The  cortex  in  a mature  spore  is  a featureless, 
electron-transparent  zone  between  tne  core  and  the  coats. 

Cortex  structure  is  related  to  type  I peptidoglycan  common  in 
vegetative  cell  walls.  Some  modifications  unique  to  spores  are 
present:  45-60%  of  the  muramic  acid  residues  lack  both  a peptide 
and  an  N-acetyl  part;  instead  they  form  an  internal  amide.  The 
core  or  protoplast  is  a relatively  normal  cell  in  terms  of  its 
raacromolecular  constituents.  Some  differences  exist  in  enzymes. 

The  properties  of  spore  ONA  are  not  significantly  different 
from  those  of  vegetative  ONA.  8acillus  subtilis  spores  contain 
a single  complete  genome,  but  other  species  may  have  multiple 
copies  (fitz-dames  and  Young,  1959).  The  physical  states  of 
enzymes  and  low  molecular  weight  solutes  in  the  core  depend  very 
critically  on  the  amount  of  water  present,  and  spores  have  lower 
water  content  than  vegetative  cells. 

The  water  activity  (Aw)  and  density  of  Bacillus  stearother- 
mophilus  spores  have  been  found  to  be  0.83  and  1.4  g/ml, 
respectively  (density  was  determined  using  a gradient  of 
tetrachloroethylene  and  ter-butyl  alcohol)  by  Algie  and  Watt 
(1984).  Other  authors  have  found  the  density  to  range  from  1.22 
to  1.34  g/ml  for  B.  stearothermophilus  and  from  1.20  to  1.35  for 
different  bands  of  8.  subtilis.  using  gradients  of  metrizamide, 
renografin  and  percoll  (Tisa  et  al.,  1982).  The  nature  of 
dehydration  during  sporogenesis  has  been  studied  by  many 


authors.  Marquis  et  al.,  (1985)  used  nondestructive  physical 
techniques  (measurement  of  the  dielectric  properties). 

Many  gram-positive  bacilli  (Including  Bacillus  and  Clos- 
tridium) form  endospores  under  certain  adverse  nutritional 
conditions  which  no  longer  allow  vegetative  growth  (Freese, 
1981).  This  process  is  called  spoliation,  and  the  cell  that 
originates  a spore  is  called  the  parent  cell.  Recently, 
sporulation  of  8.  subtilis  was  reviewed  thoroughly  by  Piggot 
(1985)  and  his  paper  is  used  here  as  the  fundamental  reference 
on  sporulation.  Exponentially  growing  (vegetative)  B.  subtilis 
can  differentiate  to  form  heat  resistant  endospores  in  response 
to  starvation  for  carbon,  nitrogen  or  phosphorus.  Differentia- 
tion of  Bacillus  subtilis  takes  approximately  six  hours  at  42 
deg.  C.  Populations  behave  fairly  synchronously,  altnough 
there  is  no  need  for  Interaction  within  the  population.  The 
basic  sequence  of  changes  during  the  formation  of  bacterial 
spores  is  similar  for  all  species  of  Bacillus  and  Clostridium 
that  have  been  examined.  The  convention  of  numbering  different 
sporulation  stages  with  Roman  numerals  is  now  in  general  use. 

The  morphological  changes  associated  with  8,  subtilis 
sporulation  are  presented  in  Figure  5.  Sporulation  can  be 
divided  into  three  successive  phases: 

a)  The  start  of  differentiation  (stages  0 to  II)  initiated 
by  the  formation  of  two  chromosomes,  and  completed  wnen  the  cell 
divides  into  two  cells  that  differ  in  size.  The  transition  from 


Figure  5.  Morphological  changes  associated 
Bacillus  subtilis  sporulation. 
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a vegetative  bacterium  to  a bacterium  including  an  asymmetri- 
cally sited  spore  septum  may  be  considered  a modified  cell 
division.  An  endogenous  factor  has  been  demonstrated  in  8^ 
cereus.  and  extracts  of  many  sporeformers  possess  the  ability  to 
transform  vegetative  cells  to  granulated  cells  having  the 
property  of  sporulating  endotropically  (Srinivasan,  1965). 

b)  In  the  second  phase  (stages  Ii  to  III),  differentiation 
becomes  fixed.  The  two  cells  formed  have  their  own  genomes  that 
may  function  differently.  When  the  septum  is  formed  at  stage 

II,  the  two  cells  differ  in  size,  but  their  contents  are 
presumably  similar.  Addition  of  rich  media  at  this  stage  causes 
the  septum  to  thicken,  and  both  cells  emerge  as  vegetative 
bacteria  (this  may  not  be  true  for  all  species).  After  stage 

III,  a mother  cell  can  no  longer  return  to  a vegetative  bac- 
terium upon  addition  of  rich  medium,  while  a forespore  can 
return  to  vegetative  growth  only  by  spore  formation  and  subse- 
quent germination.  From  stage  III  on,  the  forespore  is  no 
longer  in  contact  with  the  medium;  it  is  entirely  surrounded  by 

c)  The  characteristic  properties  of  a mature  spore  are 
developed  in  the  third  phase  (stages  III  to  VII).  Formation  of 
the  cortex  and  primordial  germ  cell  wall  between  the  opposed 
membranes  is  defined  as  stage  IV,  and  deposition  of  the  pro- 
teinaceous spore  coat  around  the  forespore  defines  stage  V.  The 
processes  overlap  in  some  species.  Synthesis  of  dipicolinic 
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acid  begins  at  about  stage  IV  in  the  mother  cell.  Calcium 
accumulates  in  forespores  slightly  prior  to  DPA  accumulation. 
Isolated  forespores  have  a specific  facilitated  diffusion  system 
for  Ca**  (Setlow,  1981).  Heat  resistance  develops  during  stage 
V and  increases  gradually.  Heat  resistance  is  characterized  by 
an  increment  in  lethal  heat  treatment  required  to  inactivate  a 
spore,  often  by  several  orders  of  magnitude  . 

Germination  is  a well  defined  stage  in  the  development  of 
spore  forming  bacteria.  It  means  conversion  of  a resistant  and 
dormant  spore  into  a sensitive,  metabolically  active  form 
(Gould,  1969).  It  is  also  described  as  a heat  labile,  non- 
refractile,  and  stainable  form,  yet  readily  distinguished  from 
vegetative  bacteria  (Keynan  and  Halvorson,  1965).  The  kinetics 
of  germination  of  aerobic  spores  were  studied  by  Vary  and 
McCormick  (1965).  The  overall  germination  process  has  been 
defined  in  terms  of  the  following  events:  a)  commitment,  b)  loss 
of  heat  resistance  and  excretion  of  calcium-DPA,  c)  loss  of 
cortex,  retractility,  and  UV  resistance  (Setlow,  1981). 

Lefrebvre  and  Antippa  (1985)  developed  a mathematical  model  of 
spore  germination  based  on  population  dynamics.  The  system  was 
defined  as  a time  ordered  sequence  of  four  distinct  phases: 
Activation,  Triggering,  Initiation,  and  Outgrowth.  Each  phase 
has  a characteristic  time  lag  and  a characteristic  rate  of 
evolution  controlled  by  a mathematical  function.  The  model 


fitted  the  experimental  data  very  closely.  This  work  is  an 
example  of  the  application  of  system  analysis  of  population 
dynamics  to  the  study  of  bacterial  spore  germination. 

Heat  Resistance  of  Bacterial  Spores 

Bacterial  spores  exhibit  the  greatest  thermoresistance  of  all 
living  organisms.  It  can  exceed  that  of  vegetative  bacteria  by 
a factor  of  105  or  more.  The  range  of  heat  resistance  between 
spores  of  different  species  is  very  large  (Gould  and  Dring, 
1974).  In  Bacillus  species  heat  resistance  generally  correlates 
with  the  maximum  growth  temperature  of  the  vegetative  form 
(Marth,  1978b).  Heat  resistance  varies  among  spore  crops,  and 
it  changes  with  composition  of  the  sporulation  medium  (Gould  and 
Bring,  1974).  It  is  not  certain  which  spore  component  or 
structure  is  the  primary  target  of  thermal  inactivation  (Gombas 
1963).  Lindsay  et  al.,  (1985)  published  a review  on  the  subject 
which  suggests  that  theories  on  heat  resistance  fall  into  two 
main  classes.  The  first  proposes  that  resistance  results  largely 
from  partial  dehydration  of  the  protoplast.  The  second  proposes 
that  heat  resistance  depends  on  molecular  rearrangements 
resulting  from  the  presence  of  calcium  dipicolinate,  glutamate, 
phosphoglycerate,  low  molecular  weight  basic  proteins  and 
possibly  other  substances  that  might  stabilize  spore  components. 
Mechanisms  discussed  in  more  detail  were  a)  reduction  in 


protoplast  volume,  b)  reordering  of  macromolecules,  c)  encase- 
ment of  the  reduced  protoplast  within  contractural ly  synthesized 
cortex  and  coat,  d)  high  concentration  of  solutes  with  con- 
comitant reduction  in  water  content  of  the  protoplast,  e) 
maintenance  of  lowered  water  content  by  pressure  developed  in 
the  cortex,  fj  development  of  a solid  support  system  of  calcium 
diplcolinate  with  specific  interactions  of  DPA/CaDPA  with 
nucleic  acids  and  possibly  proteins,  and  g)  constraint  of 
macromolecular  motions. 

Koshikawa  et  al.,  (1984)  showed  that  divestment  of  integument 
layers  exterior  to  the  cortex  had  little  influence  on  heat 
resistance  of  Bacillus  meqaterium.  Among  four  morphotypes  of 
spores  tested,  heat  resistance  did  not  correlate  with  total 
water  content,  wet  density,  refractive  index,  or  dipicolinate  or 
cation  content,  but  it  did  correlate  with  the  volume  ratio  of 
protoplast  to  protoplast  plus  cortex.  Murrell  and  Uartn  (1965) 
determined  calcium,  magnesium,  alpha-epsilon-diaminopimelic  acid 
(DAP),  dipicolinic  acid  and  hexosamine  in  20  spore  preparations 
of  8acillus  species,  with  a 700-fold  range  in  heat  resistance, 
and  examined  their  relation  to  heat  resistance.  Calcium  showed 
a small  but  significant  increase  with  heat  resistance, while 
dipicolinic  acid  did  not  Increase  significantly.  The  best 
multiple  regression  equation  for  predicting  heat  resistance 
included  DAP,  calcium  dipicolinate  and  the  Mg-Ca  ratio 


43 

(r=0.943).  Warth  (1985)  proposed  that  partial  dehydration  Is  a 
tenable,  general  mechanism  for  stabilization  of  spore  proteins 
to  heat.  An  Aw  of  0.73  was  required  to  stabilize  a spore  enzyme 
in  vitro  to  the  same  extent  found  in  vivo.  Removal  of  calcium 
dipicolinate  and  other  low  molecular  weight  spore  substances  did 
not  have  a major  effect,  thus  suggesting  that  dehydration  of  the 
core  is  maintained  by  water  tension  in  the  cortical  layers. 

Hurrell  (1981)  reported  on  biophysical  studies  of  bacterial 
spores  and  concluded  that  the  ONA-OPA-CaDPA  Interaction  and 
reduction  of  protoplast  Aw  might  be  the  main  mechanisms  of  spore 
heat  resistance.  Lindsay  and  Murrell  (1985,  1986)  studied  the 
effect  of  dipicolinic  acid  and  its  calcium  chelate  on  the 
spectral  characteristics  and  the  wet  and  dry  densities  of  DNA. 
They  concluded  that  interaction  of  dipicolinic  acid  with  ONA 
influences  the  spore  wet  density  and  the  ratio  of  core/(core- 
♦cortex)  volume,  and  it  may  influence  spore  heat  resistance. 
Dipicolinic  acid  was  found  to  displace  ethidium  bromide  from 
DNA,  indicating  that  one  mechanism  of  binding  is  intercalation 
(Lindsay  and  Murrell,  1986).  That  hydrophobic  interaction  would 
suggest  a mechanism  whereby  water  could  be  displaced  from  spore 
polynucleotides,  increasing  their  resistance  to  denaturation. 
Handley  et  al.,  (1981)  investigated  the  binding  of  ethidium 
bromide,  and  showed  strong  binding  to  spores  which  reduced  the 
heat  resistance  of  the  spores  after  exposure  to  ethidium 
bromide. 


Chandler  (1973)  studied  the  organization  of  DNA  within 
sporulating  cells  of  Bacillus  species  and  showed  that  some 
species  may  contain  more  than  one  genome  in  their  spore  (also 
Ooi,  1969).  Zytkovicz  and  Nalvorson  (1972)  reported  the 
existence  of  a dipicolinate  acid-less  mutant  of  B.  subtills 
spores  that  completely  retained  heat  resistance;  other  OPA-less 
mutants  retained  only  a fraction  or  completely  lost  heat 
resistance. 

Stewart  et  al.,  (1981)  reported  on  the  distribution  of 
calcium  and  other  elements  in  B.  cereus.  Almost  all  the 
calcium,  magnesium,  phosphorus  and  manganese  were  located  in  the 
core  region,  and  a high  concentration  of  silicon  was  found  in 
the  cortex/core  layer.  Bender  and  Harquis  (1985)  studied  spore- 
specific  mineralization.  They  concluded  that  it  is  an  important 
factor  in  spore  heat  resistance,  but  the  relationship  between 
resistance  and  mineralization  is  complex  and  dependent  on  lethal 
temperature.  Beaman  and  Gerhardt  (1986)  studied  the  relation- 
ship among  heat  resistance  and  protoplast  dehydration,  minerali- 
zation and  thermal  adaptation.  They  found  that  increases  in 
mineralization  and  thermal  adaptation  reduced  protoplast  water 
content  between  limits  of  ca.  57  and  28%  (wet  weigh  basis)  and 
correlated  with  increases  in  thermal  resistance.  Above  and 
below  those  limits,  increases  in  mineralization  and  thermal 
adaptation  correlated  with  increase  in  spore  resistance  indepen- 
dent of  unchanged  protoplast  water  content.  The  authors 


minerals  might 


suggested  a molecular  mechanism  by  which 
stabilize  spores  by  replacement  of  water  in  biopolymers,  quoting 
as  evidence  the  intercalation  of  DNA  by  calcium  dipicolinate. 

The  influence  of  water  activity  on  the  heat  resistance  of 
bacteria  was  reviewed  by  Hurrell  and  Scott  (1966),  Corry  (1975), 
Beaman  et  al.,  (1982),  Hurrell  (1981),  and  Lindsay  et 
a 1 . , ( 1985 ) . Spore  heat  resistance  as  a function  of  Aw  reaches  a 
maximum  between  0.2  and  O.S.  The  effect  of  aoded  solute  was 
also  discussed.  The  influence  of  pH,  Aw  and  recovery  tempera- 
ture on  Bacillus  spore  heat  resistances  of  151  different 
strains  was  explored  by  Braithwaite  and  Pertgo  (1971).  Between 
pH  6.4  and  9.0,  severe  heat  treatments  were  needed  to  inactivate 
spores,  outside  this  range  relatively  mild  processes  are 
sufficient.  Beaman  and  Gerhardt  (1986)  investigated  the 
correlation  of  spore  heat  resistance  with  protoplast  dehydra- 
tion, mineralization  and  thermal  adaptation.  All  three  factors 
contributed  to  and  were  necessary  for  spore  heat  resistance,  but 
dehydration  predominated.  Beaman  et  al.,  (1982)  investigated 
the  correlation  of  neat  resistance  with  water  content,  wet 
density  and  protoplast/sporoplast  volume  ratio.  An  exponential 
increase  in  heat  resistance  correlated  directly  with  wet  density 
and  inversely  with  water  content  and  protoplast/sporoplast 
volume  ratio.  Algie  (1983J  also  showed  a correlation  between 
heat  resistance  and  the  average  decrease  in  protoplasm  volume 
during  sporulation. 


resistance  of  Cl.  botulinum 


The  heat 

0.9  was  reported  by  Alderton  et  a 1 . . (1980).  The  maximum 
reached  was  ca.  Aw  = 0.4.  The  correlation  between  protoplast 
wet  density  and  protoplast  water  content  was  reported  by  Lindsay 
et  al.,  (1985)  as  a method  for  obtaining  the  protoplast  water 
content  from  the  buoyant  density  in  Metrizamide  gradients.  The 
average  protoplast  dry  density  was  1.460  g/ml.  The  protective 
effect  of  fat  on  spore  heat  resistance  was  investigated  by 
Senhaji  and  Loncin  (1977),  Senhaji  (1977),  and  Ababouch  et 
el., (1986).  The  reduction  in  Aw  was  very  important.  The  last 
authors  also  proposed  a theory  of  spore  clumping  to  explain  the 
shape  of  experimental  curves  obtained  as  one  possible  explana- 


Activation 

Activation  is  a process  that  changes  viable  spores  incapable 
of  germinating  in  an  enumeration  medium  into  activated  spores 
capable  of  germination  and  growth.  The  usual  method  to  activate 
spores  is  exposure  to  heat.  Activation  occurs  for  pH  in  the 
range  2.0  to  8.5  (Keynan  et  al.,  1965).  Activation  was  reviewed 
by  Keynan  et  al  (1965);  they  mention  that  activation  was  first 
discovered  in  mold  spores,  that  it  was  used  as  early  as  1919  by 
Weizmann,  and  that  the  main  agents  of  activation  are  exposure  to 
heat,  to  low  pH  (1-1.5),  thiol  compounds  and  strong  oxidizing 


agents.  Their  discussion  of  possible  activation  mechanisms 
supported  the  idea  that  the  macromolecules  responsible  for 
maintenance  of  the  dormant  state  are  proteins  rich  in  cystine 
and  stabilized  by  S-S  linkages.  Handley  (1977)  said  that  the 
activation  process  was  related  to  changes  in  tertiary  structure 
of  a protein  responsible  for  the  maintenance  of  dormancy. 

Freese  and  Cashel  (1965)  found  that  spores  of  8.  subtills  can 
be  activated  by  exposure  to  calcium  dipicolinate.  Levinson  and 
Hyatt  (1969)  investigated  activation  of  8.  meqaterium  spores 
and  found  that  they  can  be  activated  for  germination  on  glucose 
by  heating  in  an  aqueous  suspension  (but  not  if  heated  dry),  by 
treating  with  aqueous  ethyl  alcohol  at  30  deg.  C,  or  by  exposure 
to  water  vapor  at  room  temperature.  They  postulated  that 
activation  involves  hydration  of  a critical  spore  site,  probably 
enzymatic,  and  that  heating  and  ethyl  alcohol  modify  the 
molecular  or  physical  properties  of  liquid  water  so  as  to 
resemble  water  vapor  in  its  ability  to  reach  or  react  with  this 

Lewis  et  al.,  (1965)  published  an  important  paper  on  dormancy 
and  activation  of  bacterial  spores  in  which  they  reported  that 
the  kinetics  of  heat  activation  and  inactivation  do  not  conform 
to  a pattern  corresponding  to  two  consecutive  first-order 
reactions.  Also,  calcium-stripped  spores  lost  heat  resistance 
but  regained  it  when  replenished  with  calcium.  Spore  activation 
in  Bacillus  subt ills,  including  some  thermodynamic  properties  of 


the  system,  was  thoroughly  Investigated  by  Busta  and  Ordal 
(1964 ) . Foerster  (1983)  found  that  spores  of  thermophilic 
bacilli  could  be  activated  by  shifting  them  to  and  holding  them 
at  temperatures  below  their  optimum  growth  temperature  of  65 
deg.  C.  The  optimum  reduced  temperature  for  activation  was  30 
deg.  C.  Spores  of  B.  stearothermophilus  NG8101  activation  at  30 
deg.  C were  investigated  by  Foerster  (1985)  who  found  that 
sodium  nitrite  was  a potent  activator.  Complete  spore 
activation  was  obtained  after  si*  hours  or  less  at  30  deg.  C. 
Sodium  nitrite  was  ineffective  as  an  activator  at  65  deg.  C. 


Injured  bacterial  spores  are  characterized  by  an  inability, 
as  a result  of  subletha!  damage,  to  develop  visible  signs  of 
growth  using  standard  growth  media  and  techniques  or  even 
optimal  for  growth  by  untreated  spores.  However,  they  retain 
colony-forming  capabilities  when  mediums  or  growth  conditions 
are  modified.  Modifications  of  growth  requirements  like 
incubation  temperature,  nutrients,  germinants  and  increased 
susceptibility  to  inhibitors  are  signs  of  spore  injury  (Adams, 
1978,  Busta,  1978).  Germination  is  a multistage  process, 
possibly  with  several  pathways  for  completion  of  each  stage. 
Complete,  irreversible  Inactivation  of  one  or  more  stages  would 
trap  a spore  protoplast  within  the  structures  maintaining 


dormancy,  and  the  spore  effectively  would  be  dead.  If  this 
inactivation  can  be  by-passed,  the  spore  is  classified  as 
Injured. 

Hany  of  the  manifestations  of  injury  fall  into  four  major 
classes:  a)  requirements  by  survivors  for  nonnutrient  germina- 
tion stimulants,  b)  modified  optimum  incubation  temperature,  c) 
Increased  sensitivity  of  survivors  to  inhibitors,  and  d)  altered 
nutritional  requirements  of  survivors  (Adams,  1978).  A general 
review  of  documented  cases  of  bacterial  spore  Injury  was 
published  by  Foegeding  and  Busta  (1981).  Heat,  irradiation, 
ethylene  oulde,  cold,  and  other  treatments  were  discussed 
together  with  the  mechanisms:  a)  membrane  damage,  b)  QNA  damage, 
c)  damaged  or  altered  germination  systems,  and  d)  damage  to 
dehydration  mechanism.  A review  of  factors  important  to  the 
detection  of  injured  bacterial  spores  in  processed  foods  was 
published  by  Johnson  and  Busta  (1984). 

The  case  of  Bacillus  subtilis  spore  Injury  has  been  investi- 
gated by  several  authors.  Edwards  et  al.,  (1965b)  studied 
injury  at  ultrahigh  temperatures.  Spores  subjected  to  ultrahigh 
temperature  may  be  Injured  but  not  completely  inactivated; 
consequently,  heat  Injury  may  be  an  Important  factor  in  the 
evaluation  of  a successful  sterilization  process.  Three  criteria 
indicated  Injury:  a)  significant  reduction  in  survivors  when 
spores  were  enumerated  with  a standard  medium,  but  not  when  the 
standard  medium  contained  added  CaCly  and  sodium  dipicolinate. 
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b)  more  survivors  of  a damaging  heat  treatment  were  enumerated 
with  the  standard  medium  after  incubation  at  32  deg.  C than  at 
45  deg.  C (opposite  to  results  obtained  with  untreated  or 
slightly  heated  spores),  and  c)  apparent  number  of  survivors 
increased  during  an  initial  period  of  storage  at  3 deg.  C when 
enumerated  with  the  standard  medium  at  45  deg.  C,  indicating 
reversible  spore  injury.  Adams  and  Busts  (1972)  presented 
results  on  heat  injury  in  B.  subtills  as  the  inactivation  of  a 
spore  germination  system.  The  thermodynamic  properties  deter- 
mined suggested  that  the  mechanism  of  injury  in  this  case  was 
protein  denaturation. 

Gould  (1984)  and  Hurst  (1984)  reviewed  repair  in  terms  of 
existing  knowledge  on  injury.  Paterson  and  Gentner  (1984) 
analyzed  mechanisms  of  ONA  repair.  The  two  principal  pathways 
for  repairing  damaged  DNA  were  a)  direct  reversal  of  an  altered 
component  back  to  its  normal  structure,  and  b)  excision  of 
defective  material  along  with  adjoining  nucleotides  followed  by 
insertion  of  a normal  nucleotide  sequence  utilizing  an  undamaged 
complementary  DNA  chain  for  base-pairing  instructions.  Heat 
induced  damage  to  bacterial  chromosome  and  its  repair  was  the 
subject  of  a review  by  Pellon  and  Sinskey  (1984). 

Oelayed  germination  of  heat  damaged  Cl.  botulinum  spores  was 
studied  by  Lynt  et  al.,  (1983).  In  that  paper,  any  culture 
requiring  more  than  one  week  to  show  evidence  of  growth  was 
considered  delayed.  Incubation  was  prolonged  up  to  13  weeks. 
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The  value  of  decimal  reduction  time  D increased  as  much  as 
fourfold,  and  proteolytic  strains  showed  suOstantial  increases 
in  2 values.  The  recovery  of  Bacillus  stearothermoohilus  used 
as  a bioindicator  organism  was  investigated  by  Caputo  et  al., 
(1980).  Three  factors  were  shown  to  affect  recovery  of  the 
biological  indicators:  a)  Incubation  temperature,  b)  post- 
sterilization hold  time,  and  c)  post-  sterilization  storage 
temperature.  These  factors  should  be  adequately  controlled  in 
order  to  obtain  reproducible  results. 

System  Analysis  of  Population  Dynamics 

System  analysis  is  a body  of  concepts  and  methods  for 
analyzing  complex  entities.  Its  development  and  application 
continually  expand.  It  is  an  attempt  to  establish  knowledge  and 
understanding,  and  its  application  to  the  fields  of  food 
engineering  and  microbiology  will  facilitate  the  solution  of 
complex,  practical  problems.  The  last  section  of  this  review  is 
a brief  examination  of  the  most  important  concepts  of  system 
analysis  pertinent  to  biological  systems.  The  main  reference 
used  for  this  section  is  Smerage  (1985). 

Because  of  the  interdisciplinary  character  of  system  analysis 
there  are  necessarily  numerous  definitions.  A system  is  a 
composition  of  objects  which  comprise  a cohesive,  functional 
whole  exhibiting  characteristic  behavioral  properties.  To  be  a 
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composition,  a system  must  consist  of  a)  a set  of  objects  and  bj 
a set  of  interconnections.  External  stimuli  received  by  a 
system  are  called  inputs,  and  the  responses  of  the  system  are 
called  outputs.  Within  a system  certain  attributes,  called 
state  variables,  measure  the  internal  condition  of  or  state  of 
the  system.  Behavioral  analysis  examines  causal  relationships 
between  input,  state  and  output  variables  to  determine  and 
characterize  behavioral  properties  of  a system. 

Populations  and  physiologies  are  amenable  to  the  same  repre- 
sentational theory  as  physical  systems  when  analyzed  from  the 
process-oriented  approach  followed  here.  Four  genera  of  process 
—storage,  transport,  transformation  and  source  or  sink— and 
their  generic  descriptions  are  adequate  to  describe  most,  if 
not  all,  types  of  systems. 

Diagrams  are  commonly  used  to  represent  and  to  support  verbal 
descriptions  of  systems.  Only  three  elements  of  the  graphical 
system  employed  in  this  research  need  to  be  defined  here;  those 
elements  are  depicted  in  Figure  6.  A store,  depicted  by  the 
object  in  Figure  6a,  represents  a population.  The  object  in 
Figure  6b  represents  a transformation  such  as  metamorphosis  or 
mortality,  of  Individuals  between  specified  populations.  A sink 
is  a population  that  is  no  longer  able  to  participate  in  the 
process  of  a system;  it  is  depicted  by  the  object  in  Figure  6c. 

Each  object  in  Figure  6 has  one  or  more  extensions  called 
terminals.  Terminals  provide  the  graphical  means  for 


Figure  6.  Representation  of  system  components: 

store  (a),  transformation  (b),  and  sink  (c). 


1 


N 
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interconnecting  component  processes  composing  a system  and  for 
representing  resulting  interactions.  Associated  with  each 
terminal  are  conjugate  variables  denoting  extensive  and 
intensive  attributes  relevant  to  the  description  of  the  object. 
The  process-oriented  approach  to  modeling  employed  here 
identifies  the  extensive 

variables  of  a terminal  with  the  flow  of  individuals  incident  to 
a component  conjugate  to  the  intensive  property  associated  with 
the  terminal. 

A population  is  a distinct  group  of  organisms  with  common 
type.  Typically  the  size  of  a population  fluctuates  with  time 
in  response  to  natality,  development,  mortality  and  migration. 
Populational  considerations  do  not  include  details  about 
individual  members,  but  attention  focuses  on  properties  of  the 
aggregate.  Extensive  properties,  such  as  number  of  individuals, 
space  occupied,  and  flows  of  individuals  depend  on  tne  size,  or 
extent,  and  content  of  a system  component.  When  two  stores  of  a 
species  are  aggregated,  extensive  parameters  of  the  composite 
store  equal  sums  of  corresponding  parameters  of  the  individual 
stores.  Intensive  properties,  such  as  temperature,  concentra- 
tion and  population  density,  do  not  depend  on  the  size  or  extent 
of  tne  component.  A principal  population  sink  is  that  for  the 
dead,  and  in  this  work,  inactivated  spores  are  stored  in  this 


In  general,  the  diagrams,  called  process  diagrams,  indicate 
component  processes,  their  manner  of  interconnection  and 
Important  system  variables.  Junctions  of  component  terminals 
Identify  interfaces  between  the  components  involved.  Interfaces 
are  sites  of  component  interactions,  and  they  are  indicated  in 
diagrams  by  dots  at  junctions,  which  are  called  vertices,  nodes, 
or  interfaces. 

There  are  several  examples  of  biological  transformations 
metamorphosis  of  insects  and  transformations  between  distinct 
living  states,  activation,  inactivation  and  injury  of  bacterial 
spores,  etc.  Kortality  transforms  individuals  from  a live  state 
to  a dead  state  regardless  of  the  cause  (physiological  failure, 
disease,  predation,  parasitism  etc.).  Development  of  an 
organism  is  a complex  physiological  and  morphogenetic  process 
where  age  classes  may  be  considered,  together  with  the  physiolo- 
gical rates  of  development  (aging).  Reproduction  is  the  key 
to  propagation  and  maintenance  of  a species.  Reproduction  is  a 
transformation  that  connects  tne  adult  population  with  the  first 
age-class  population.  Fecundity  implicitly  depends  on  environ- 
mental factors,  and  that  dependence  may  be  accounted  for  adding 
related  terminals  to  the  transformation  object.  A transforma- 
tion representing  mortality  removes  members  from  a store  of  live 
individuals  and  transfers  them  to  the  populational  sink. 
Predation  and  parasitism  are  distinct  forms  of  mortality,  but 
their  representation  and  discussion  is  similar:  a transformation 


representing  predation  (parasitism)  expresses  the  causal 
dependence  of  prey  (host)  flow  on  prey  (host)  density  and 
predator  (parasite)  density;  terminals  of  the  transformation 
connect  the  affected  store  of  prey  (host)  species  to  a sink  for 
the  dead  and  implicate  the  store  of  relevant  predatory  (parasi- 
tic) species. 

Hodels  are  operational  statemencs  of  knowledge  and  their 
application  strengthens  different  fields  of  knowledge.  Hodels 
of  populational  systems  are  developed  by  logically  and  systema- 
tically performing  four  steps. 

1.  Develop  a conceptual  model  by  decomposing  a given  real 
system  into  a network  of  perceived  populational 
components. 

2.  Experimentally  and/or  theoretically  determine  the 
description  of  each  component  in  the  conceptual  model. 

3.  Formulate  the  structural  description  of  the  system,  as 
depicted  in  the  conceptual  model. 

4.  Derive  a mathematical  model  corresponding  to  tne 
conceptual  model  and  descriptions  of  components  and 
structures.  It  will  consist  of  one  or  more  equations 
representative  of  the  whole  system  as  an  entity. 

Generic  descriptions  of  the  form  fundamental  families  of 
processes  guide  experimental  and  theoretical  descriptions  of 
specific  real  processes.  The  laws  of  system  structure  facilitate 
mathematical  description  of  that  aspect  of  a system  (Smerage, 


53 

1979).  Several  methods,  including  state  variaoles,  circuit  and 
node  formulations,  e»1st  for  formulating  a mathematical  model  of 
a system  (Koenig  et  al.,  1967,  Martens  and  Allen,  1969). 


MATERIALS  AND  METHODS 


The  scope  of  work  undertaken  in  this  study  has  been  divided 
into  two  parts  addressing  the  two  most  important  factors 
■affecting  validation  of  commercial  sterilization  processes  with 
bioindicator  units  (Bills).  These  factors  are  the  heat  transfer 
characteristics  of  the  unit  and  the  Inactivation  kinetics  of  the 
indicator  microorganism. 

A series  of  experiments  was  performed  to  determine  the 
contributions  of  heat  transfer  and  material  strength 
characteristics  of  Bills  to  the  discrepancy  between  predictions 
by  3igelow's  kinetic  model  and  corresponding  experimental  data. 
In  these  experiments,  survivor  3.  stearothermoohilus  spores  from 
different  heat  treatments  were  enumerated  and  conpared  to 
corresponding  expected  values  calculatea  using  traditional 
kinetic  models  with  the  time-temperature  history  of  the 
treatments  . Later,  Bills  were  designed  and  built  from  aluminum 
for  improved  heat  transfer  and  strength  characteristics,  and 
their  performance  was  compared  to  that  of  plastic  Bills  normal  ly 
used  in  commercial  practice  . 

Secondly,  the  kinetics  of  bacterial  spore  inactivation  were 
studied  in  four  steps.  In  step  one,  a new  kinetic  model  was 
developed  of  a suspension  of  bacterial  spores  treated  at  lethal 
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temperatures.  The  conceptual  model  was  proposed  first  from  an 
application  of  tne  theory  of  population  dynamics,  which  led  to 
development  of  a mathematical  model  as  a set  of  differential 
equations  describing  the  various  reactions  (transformations) 
involved.  The  solution  of  these  differential  equations  produced 
response  equations  that  consisted  of  sums  of  exponential  terms. 

The  next  step  in  this  part  of  the  work  was  determination  of 
parameters  of  the  model;  it  consisted  of  laboratory  experiments 
involving  carefully  controlled  time-temperature  exposures  of  8. 
subtilis  spore  suspensions.  Data  from  these  experiments  were 
used  to  produce  experimental ly-derived  survivor  curves  for  both 
normal  and  injured  spores.  These  curves  were  analyzed  to  deter- 
mine values  of  parameters  of  the  model,  which  were  coefficients 
and  exponents  of  exponential  terms  in  the  response  equations. 
Subsequently,  predicted  survivor  curves  were  compared  witn 
experimental  survivor  curves  at  each  lethal  temperature  used  in 
the  thermal  exposure  experiments. 

In  the  third  step,  temperature  dependence  of  the  first  order 
reaction  rate  constants  was  determined  in  accordance  with  the 
Arrhenius  equation,  and  it  substituted  for  the  rate  constants  in 
the  new  kinetic  model.  In  this  form,  the  model  should  be 
capable  of  predicting  the  number  of  survivors  of  a thermal 
process  in  which  the  lethal  temperature  was  varied  as  an 
arbitrary  function  of  time  (a  transient  temperature  process). 


The  fourth  step  was  a final  test  or  validation  of  the  model 
in  which  experimental  results  of  a transient  temperature  process 
were  compared  with  results  predicted  by  the  model  for  the  same 
temperature  history. 

Oetails  of  methodology  are  presented  in  the  following 
subsections,  beginning  with  a description  of  how  spore  suspen- 
sions were  prepared  and  handled  throughout  tne  study.  Then, 
the  methodology  used  in  carrying  out  the  heat  transfer 
experiments  with  bioindicator  units  is  discussed.  Finally, 
development  and  testing  of  the  new  kinetic  model  are  discussed. 

Theory  formation  provided  the  rationale  for  the  methodology 
described  in  the  following  sections,  identifying  significant 
variables,  tneir  ranges  of  variation  and  tne  corresponding 
procedural  aspects.  Hence,  theory  formation  led  to  a logically 
structured  experimental  design.  It  is  presented  for  the  study 
of  both  heat  transfer  in  BIUs  and  kinetics  of  bacterial  spore 
inactivation  in  terms  of  available  Information  gatnered  from  the 
literature  review  to  support  various  assumptions  and  expected, 
logical  inferences. 

Spore  Suspensions  and  Enumeration  of  Survivors 


Bacillus  stearothermophilus  spores.  The  spore  suspension 
used  was  provided  by  Ross  laboratories,  already  placed  inside 
Bills.  The  thermal  resistance  parameters  were  independently 
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determined  as  O250  s 6.73  min  and  2 = 16.16  deg.  F by  an 
outside  laboratory  at  the  University  of  Massachusetts  (Evans, 
1984).  The  suspension  was  kept  at  4 deg.  C. 

Bacillus  subtilis,  strain  "A"  spores.  The  suspension  used 
was  more  than  three  years  old,  and  was  obtained  from  Prof.  F.F. 
Busta  at  the  University  of  Florida  Food  Science  and  Human 
Nutrition  Department,  washed  and  stored  at  4 deg.  C.  Inspection 
under  the  phase  contrast  microscope  showed  the  spores  to  be 
fully  refractile  and  free  from  vegetative  cells,  phase  dark 
spores  and  clumps.  The  direct  microscopic  count  was  performed  in 
a phase  contrast  microscope  using  a Bright-Line  hematocytooeter 
(American  Optical). 

Enumeration  of  Survivor  Soores.  For  the  enumeration  of  B. 
stearothermophilus  survivor  spores,  the  growth  medium  was 
Spectrum  soy  agar,  provided  by  Ross  laboratories.  Tne  medium 
was  prepared  by  mixing  20  g in  500  ml  of  distilled  water,  and  pH 
was  adjusted  to  7.3.  Then  the  medium  was  placed  in  free  flowing 
steam  for  30  min  and  stirred  and  sterilized  at  250  deg.  F (121.1 
deg.  C)  for  20  min.  Oecimal  dilutions  of  the  experimental 
samples  were  prepared  in  sterile  distilled  water,  and  1 ml  of 
each  dilution  was  placed  in  a Petri  plate  with  c°.  25  ml  of 
growth  medium  added.  The  incubation  temperature  was  55  deg.  C 
for  48  hr. 

The  growth  medium  used  for  B,  subtilis  was  that  described  by 
the  case  when  no  calcium  dipicolinate 


Edwards  et  al.,( 1965a),  ir 
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»as  added.  Two  incubation  temperatures  were  used:  45  deg.  C for 
18  hrs.  (standard  incubation  temperature  for  tnis  strain),  and 
32  deg.  C for  72  hrs.  (for  detection  of  injured  spores).  Study 
of  the  type  of  spore  injury  measured  in  this  work  is  dependent 
upon  simultaneous  or  independent  activation  of  a germination 
mechanism  tnat  allows  spores  to  grow  at  lower  temperatures 
(Edwards.  1965). 

A test  tube  mixer  (Vortex-Genie)  was  used  to  stir  dilution 
tubes  and  BIUs  for  15  s.  In  all  cases,  decimal  dilutions  of  the 
experimental  samples  were  prepared  in  sterile  distilled  water, 
and  1 ml  of  each  dilution  was  placed  in  a Petri  plate  with  ca. 

25  ml  of  growth  medium  added.  In  the  case  of  B.  subtills,  two 
plates  were  prepared  of  each  dilution,  corresponding  to  the  two 
incubation  temperatures  used.  When  the  medium  solidified  the 
plates  were  incubated,  and  the  resulting  colonies  were  counted. 

Bioindicator  Units 


In  this  section,  procedures  corresponding  to  the  study  of 
BIUs  are  presented.  Propositions,  assumptions  and  logical 
inference  that  led  to  the  experimental  design  are  presented 


Propositions 

1.  The  temperature  history  of  a bioindicator  unit  should  be 
as  close  as  possible  to  that  of  the  cold  spot  of  the  product 
(Joslyn,  1983,  Bruch,  1974). 

2.  Pressure  changes  can  affect  the  lethality  (Teiteira  et 
al.,  1985,  Cleland  and  6esterltamp  1983). 

3.  The  time-temperature  history  can  be  measured  using 
copper-constantan  thermocouples. 

4.  The  number  of  survivor  spores  can  be  measured  microbiolo- 
gically  (Pflug  et  al.,  1980). 


1.  The  conductivity  of  the  material  used  to  fabricate  a BID 
plays  an  important  role  in  its  ability  to  accurately  transmit 
the  temperature  history  of  the  product  undergoing  thermal 
treatment. 

2.  The  mechanical  strength  of  a bioindicator  unit  is 
important  when  pressure  changes  occur  during  a sterilization 

3.  The  effect  of  container  material  on  the  performance  of  a 
bioindicator  unit  can  be  studied  using  heat-penetration 
measurements  and  microbiological  enumeration  of  survivor  spores. 
Inference 

Improved  performance  of  bioindicator  units  can  be  achieved 
when  fabricated  from  materials  with  high  thermal  conductivity 


and  mechanical  strength  (i.e.  metal  vs.  plastic)  and  this 
improvement  is  more  important  for  rapid  convection-heated  liquid 
products  than  for  slow  conduction-heated  products. 

Experimental  Work 

Following  the  guidelines  established  in  theory  formation, 
plastic  BIUs  were  tested  using  both  constant  and  variable 
temperature  processes.  Next,  plastic  BIUs  were  compared  to 
aluminum  BIUs  in  a constant  temperature  process.  The 
experiments  with  plastic  BIUs  were  designed  to  test  tne 
traditional  kinetic  model  under  different  process  conditions, 
while  the  comparative  experiments  using  plastic  and  aluminum 
BIUs  at  constant  retort  temperature  were  designed  to  test  the 
relative  importance  of  heat  transfer  to  the  valioation  of 
thermal  sterilization  processes  using  81  Us . The  higher  thermal 
conductivity  of  aluminum  snould  allow  the  internal  temperature 
of  an  aluminum  BIU  to  follow  the  product  temperature  more 
closely  and,  thereby  minimizing  inaccuracies  related  to  heat 
transfer. 

Bioindicator  units.  The  plastic  BIUs  (Pflug  et  al.,  1980) 
were  provided  by  Ross  Laboratories  (Columbus,  OH).  The  aluminum 
units  were  made  with  the  same  dimensions  as  the  plastic  BIUs, 
using  a two-part  construction  in  which  the  end  section  holding 
the  spore  suspension  was  made  of  aluminum,  while  the  remaining 


66 

support  section  was  made  of  plastic  to  keep  the  aluminum  section 
well  insulated  from  tlie  can  wall-steam  interfdce.  Some  Bids  were 
Instrumented  with  needle  thermocouples  and  contained  distilled 
water.  A photograph  of  the  BIUs  is  presented  in  Figure  7. 

Procedure.  The  heat  treatments  for  experiments  with  plastic 
and  aluminum  Bills  were  performed  using  a vertical  stii  i cook 
retort.  Figure  8 shows  a 46-ounce  can  (0.1778x0.1048  m) 
equipped  with  copper-constantan  needle  thermocouples  (O.F. 
Ecklund  custom  thermocouples.  Cape  Coral,  Fla.]  and  BIUs.  The 
different  probes  were  placed  at  approximately  the  same  position 
inside  the  can.  The  time-temperature  data  were  recorded  on  a 
Digi strip  II  data  logger  (Kaye  Instruments,  Bedford,  Hass.) 
connected  to  a RT-11  computer  system  (Digital  equipment 
corporation,  Haynard,  Hass.].  The  cans  prepared  with  BIUs  and 
thermocouples  were  filled  with  hot  distilled  water  until  there 
was  a 0.2  in.  headspace;  tnen  they  were  sealed  and  placed  inside 
the  retort.  The  thermocouples  were  connected  to  the  temperature 
recorder,  and  the  retort  was  closed.  After  the  selected  venting 
time  of  2 min,  the  venting  valve  was  closed,  and  the  temperature 
was  regulated  to  follow  the  prescribed  thermal  history.  When 
the  process  time  was  reached,  the  steam  valve  was  closed;  the 
pressurieeo  air  valve  was  opened,  and  cooling  water  was  allowed 
to  enter  the  retort.  Afterwards,  the  retort  was  opened,  and  the 
3IUs  were  recovered  from  the  cans.  The  BIUs  were  immediately 
placed  in  ice  and  stored  under  refrigeration. 


Fiaure  7.  Plastic  and  aluminum  bioindicator 
’ units,  one  of  them  instrumented  with 

a needle  thermocouple. 
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Figure  8. 


Can  (46  ounces)  equipped  with  a copper 
constantan  nr-4’"  ,h“''mnrmlDle  and 
plastic  ai 


Lethality  was  calculated  by  numerical  integration  of  the 


following  standard  expression: 


t (T-250J/Z 

F = J « (14) 

where  F is  the  lethality,  T is  the  temperature  at  each  time 
interval  sampled  in  deg.  F,  t is  time  in  min,  and  Z is  the 
temperature  dependency  factor  for  thermal  resistance 
(temperature  increment  for  a ten-fold  change  in  lethal  rate. 
Temperature  history  data  from  the  heat  treatments  were  used  in 
these  integrations.  The  calculated  lethality  was  compared  with 
the  lethality  found  by  microbiological  enumeration  of  the 
survivor  spores,  which  can  be  expressed  as  follows: 

F = D250  log(No/M) (15) 

where  N and  No  are  final  and  initial  numbers  of  survivor  spores 
respectively  and  O250  is  the  decimal  reduction  time  for  the 
spore  suspension  at  the  reference  temperature. 

During  tne  course  of  this  work,  it  became  apparent  that 
differences  in  mechanical  strength  between  the  plastic  and 
aluminum  might  be  a factor  in  explaining  observed  behavior. 
Materials  testing  was  performed  to  study  high  temperature 
softening  of  plastic  Bids.  Force-deformation  curves  for  the 
plastic  Bids  at  high  temperatures  were  determined  using  a MTS 


Material  Test  System  (Houston  Instruments)  equipped  with  an 
omn i graph i c 2000  recorder.  Limitations  of  the  environmentally- 
controlled  chamber  did  not  permit  experiments  at  higher 
temperatures.  The  samples  were  hollow  cylinders  prepared  from 
the  plastic  Bills  with  0.253  in.  (0.00643  m)  outside  diameter, 
0.125  in.  (0.003175  m)  internal  diameter,  and  0.255  in. 

(0.006477  m)  height.  The  speed  was  0.0833  mm/s.  Two  tests  were 
performed:  one  at  70  deg.  C and  another  at  93  deg.  C. 

The  plastic  Bills  were  subjected  to  two  thermal  regimes:  a 
constant  retort-temperature  of  250  deg.  F (121.1  deg.  C)  and  a 
variable  temperature  wherein  the  retort  temperature  was  raised 
first  to  250  deg.  F (121.1  deg.  C),  and  then  to  260  deg.  F 
(126.67  deg.  C).  The  type  of  process  was  initially  considered 
as  an  important  factor  in  explaining  the  divergences  between 
lethality  found  by  integration  using  Bigelow's  model  and 
lethality  found  by  enumeration  of  survivors.  The  experiments  to 
compare  the  performance  of  aluminum  with  plastic  BIUs  were 
performed  at  a constant  temperature  of  250  deg  F (121.1  deg.  C). 
Therefore,  the  first  set  of  experiments  attempted  to  test  the 
effect  of  different  types  of  heat-treatment  on  the  performance 
of  plastic  Bills,  while  the  second  set  tested  the  importance  of 
the  heat  conductivity  and  mechanical  strengtn  on  the  performance 


different  BIUs. 


spore  inactivation 


In  this  section,  the  propositions,  assumptions  and  logical 
inference  in  the  theory  formation  and  the  systematic  development 
of  a model  of  a suspension  of  bacterial  spores  are  presented. 

The  experiments  performed  for  parameter  determination  and 


validation  of  t 

:he  model  are  described. 

Propositions 

1.  Activation,  inactivation,  and  injury  are  concurrent 
reactions  that  occur  during  heat  treatment  of  bacterial  spores 
(Busta,  1978;  8usta  and  Ordal,  1964;  Edwards,  ml.,  1965a). 

2.  Systems  analysis  is  a set  of  techniques  that  can  be 
applied  to  the  study  of  population  dynamics  with  multiple 
concurrent  reactions  (Smerage,  1985). 

3.  The  number  of  survivor  spores  can  be  measured  microbiolo* 
gically  (Stumbo,  1965). 

4.  The  temperature  history  of  heat  treatments  can  be 


precisely  measu 

ired  and  controlled  (Oatta  et  al.,  1986). 

Assumptions 

The  transformations  can  be  described  by  pseudo  first- 


kinetic  models. 


2.  The  transformations  can  occur  simultaneously. 

3.  The  transformations  are  independent  of  each  other. 

4.  The  dependence  of  each  rate  constant  on  temperature  can 
be  described  by  an  Arrhenius-like  moael. 

Inference 

The  effect  of  lethal  temperatures  on  bacterial  spores  can  be 
modeled  using  systems  analysis  of  population  dynamics,  if  the 
microbiological  enumeration  of  the  survivor  spores  and  the 
measurement  of  the  time-temperature  history  can  be  accurately 
determined  for  a given  spore  suspension. 

Description  of  Research  Work 

In  this  phase  of  work , and  in  accordance  with  theory 
formation,  the  kinetics  of  spore  survival  at  lethal  temperatures 
were  studied  as  a major  factor  in  the  biological  validation  of 
thermal  sterilization  processes.  The  main  changes  (transforma- 
tions) that  occur  during  exposure  were  identified  from  a review 
of  the  literature  and  incorporated  in  a conceptual  model 
amenable  to  systematic  analysis.  A mathematical  model  was 
generated  and  a series  of  experiments  at  different  constant 
temperatures  were  performed.  These  experiments  provided  data 
for  estimation  of  parameters  in  the  model,  including  Arrhenius 


description  of  the  temperature  dependence  of  rate  constants.  As 
an  initial  validation  of  the  model,  the  results  of  a 
nonisothermal  experiment  were  compared  with  predictions  of  the 
model  for  the  temperature  history  recorded  during  the 
experiment. 

Model  Development 

A conceptual  model  of  a suspension  of  bacterial  spores  at 
lethal  temperatures  was  proposed,  and  components  and 
interactions  were  identified.  The  components  of  the  system  are 
the  different  stages  of  a population  of  bacterial  spores 
undergoing  a heat  treatment  at  lethal  temperatures:  mature 
nonactivated  spores,  activated  spores,  injured  spores,  and  a 
less  heat  resistant  population.  The  transformations  under 
consideration  are  inactivation  (0),  activation  (A)  and  injury 
(I).  Stores,  representing  the  populations,  and  transformations 
were  analyzed  together  with  their  possible  interactions. 
Mathematical  expressions  describing  the  transformations  were 
then  incorporated,  and  a system  of  differential  equations  in  the 
subpopulations  was  generated.  Solution  of  this  system  of 
differential  equations  for  the  isothermal  case,  taking  into 
consideration  the  initial  conditions,  allowed  prediction  of  the 
populations  at  different  moments. 
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Conceptual  model.  The  measurement  procedure  is  very 
important  to  the  definition  of  components  in  the  system:  Spores 
are  considered  active  if  they  are  able  to  grow  when  incubated  at 
45  deg.  C.  They  are  considered  injured  if  they  are  capable  of 
growth  when  Incubated  at  32  deg.  C,  but  not  at  45  deg.  C.  They 
are  considered  inactivated  if  they  no  longer  are  able  to 
germinate  and  grow  when  incubated  at  eitner  temperature.  A less 
heat  resistant  population  consists  of  germinated  spores  that 
have  lost  their  heat-resistance  or  vegetative  cells  or,  in  the 
case  of  multiple  genome  bacteria,  spores  that  are  genetically 
less  heat-resistant.  However,  a less  heat-resistant  population 
can  produce  colonies  in  an  enumeration  medium  that  are 
indistinguishable  from  colonies  originated  by  mature  activated 
spores  of  the  system. 

The  conceptual  model  is  presented  diagrammatical iy  in  Figure 
9.  Mature  spores  that  have  not  been  activated  (dormant)  are 
potentially  able  to  produce  colonies  after  activation;  they 
comprise  the  population  in  store  Nj.  Activated  spores,  immeaia- 
tely  capable  of  producing  bacterial  colonies  in  the  enumeration 
medium,  are  represented  by  store  Ng.  Injured  spores  are 
determined  by  their  ability  to  grow  at  relatively  low  incubation 
temperature  (32  deg.  C);  tney  are  denoted  by  tne  additional 
subscript  i.  An  injured  population  has  analogous  subpopulations 
of  mature  nonactivated  and  activated  spores:  Njj,  and  Njy 
respectively,  in  this  case,  the  germination  mechanism  allows 


Figure  9.  System  diagram  for  a oacterial  s 
suspension  under  treatment  at  le' 
temperatures.  The  subscript  3 


subscript  I 


resistant  fraction,  the 
refers  to  viable  dormant 
subscript  2 refers  to 

subscript  i 
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mechanism  or  injury. 


activated  spores  to  grow  at  lower  incubation  temperatures.  The 


populations  of  less  heat-resistant  microorganisms  that  may  be 
present  are  represented  by  stores  N3,  and  Nj3  for  normal  and 
Injured  classes  respectively.  Another  component  of  possible 
Importance  , but  not  included  in  this  work,  is  the  fraction  of 
spores  that  may  mature  during  the  test  Interval,  thereby 
increasing  the  number  of  spores  in  Nj  or  iij.  This 
transformation  would  lead  to  dispersion  of  the  values  obtained 
for  the  number  of  survivors.  It  should  be  minimized  by  using  a 
homogeneous  spore  suspension. 

The  hypotheses  that  the  transformations  are  simultaneous  and 
independent  were  very  important  to  the  conceptual  model.  It  has 
four  major  independent  parts.  The  hypothesis  of  consecutive 
reactions  has  been  previously  studied  and  found  incorrect  (Lewis 
et  al.,  1965),  and  no  reason  could  be  found  to  justify  further 
coupling  of  the  reactions.  The  upper  left  part  of  Figure  9 
represents  inactivation  of  a less  heat-resistant  fraction  H3. 
The  upper  right  part  represents  activation  («2)  of  mature  viable 
dormant  spores  (Nj)  and  inactivation  of  both  populations  in 
agreement  with  the  hypotheses  that  these  transformations  are 
simultaneous  and  independent.  The  lower  part  of  the  diagram 
represents  injury-related  events  in  an  analogous  manner:  the 
left  part  represents  the  inactivation  of  a less  heat-resistant 
population  (»j3)  that  could  be  enumerated  using  a lower 
incubation  temperature  (32  deg.C);  the  right  part  represents  the 


Injury  of  viable  spores,  together  with  the  inactivation  of  both 
populations  (Nn,Ni2-  This  type  of  injury  may  be  described  as 
the  activation  of  a germination  mechanism  that  operates  at  lower 
temperatures  than  the  normal  germination  mechanism.  The 
triangles  represent  the  sink  for  the  inactivated  spores. 

Mathematical  model.  Since  the  four  parts  of  the  model  are 
not  coupled  they  can  be  analyzed  independently.  The  components 
are  described  as  follows: 


.(16) 


•(17) 


Transformations: 

Activation: 

12  * KaNi;  (18) 

Inactivation: 


Tjd'MMj;  (J  = 1.  Z)  (19) 

Inactivation  of  less  heat-resistant  fraction: 


(20) 


'3d  = <3oN31 


flows  between 


mathematically 


on  the  law  of  conservation  of  matter:  since  the  node  can  not  act 
as  a store,  the  summation  of  flows  at  every  node  is  aero. 
Consequently,  the  structure  of  the  system  is  described  by  the 
following  equations: 

Summation  of  flows  at  node  a: 

- \ ' Yld  - Y12  ' 0 or,  Tl  ■ - ( Yld  * Y12)  •••(2D 
Summation  of  flows  at  node  b: 

Id  - V2  - \ * 0 or:  Y2  =■  YJ2  - Y2q (22) 

Substitution  of  the  equations  for  the  transformations 
(18,  19,  20)  in  Equations  21  and  22  yields 

Y1  = - (KdNi  ♦ KaN!)  * - (Kd  * Ka)  Nj  (23) 

Y2  = Kali!  - KdN2  (24) 

The  mathematical  model  corresponding  to  this  part  of  the 
conceptual  model  is  obtained  by  substituting  Equations  23  and  24 
in  the  equations  for  the  stores  (16  and  17).  It  is  presented 
with  the  initial  conditions  used  in  this  research: 
dNt 

» - (Kd  * Ka)  Mi;  Ni  (0)  = Nig  ...(25) 


KaN]  - KdN2; 


N2  (0)  = N20  ...(26) 
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Solution  of  the  model  with  the  initial  conditions  yields 
the  response  equations: 

Ni  (t)  - Hi0e*p(-(Kd  * Ka)t)  (27) 

H2  (t)  = (Nio-N2o)e*p(-Kd  t)-Nioexp(-(Kd  * Ka)t)  ..(28) 

By  similar  means,  a single  differential  equation  is  derived 
for  the  less  heat-resistant  fraction,  together  with  an  initial 
condition: 

— * Kd3N3:  li3  (0)  = H3o  (29) 

Tne  solution  is 

N3  = N3oexp(-Kd3  t)  (30) 

The  case  of  injured  spores  or  activation  of  an  alternate 
germination  mechanism  was  handled  similarly.  The  mathematical 
model  is  presented  in  Table  1,  and  the  corresponding  equations 
for  the  responses  are  presented  in  Table  2. 


pension  of  bacterial 


<ia3Hi3 


"13  (0) 


Table  2.  Response  equations  of  the  model  for  constant 
temperature. 

Nj  (t)  = HioexpHKd  * Ka)t) 

N2  (tj  = (N10  ♦ Njo)exp(-Kd  t)  - Ni0exp(-(Kd  ♦ Kajt) 

N3  (t)  ■ N3oexp(-Kd3  t) 

Nil  (t)  ■ Nj10exp(-(K1d  * KiaJtJ 

Niz  ct)  ■ (Njio  * Hi2o)exp(-Kid  t)-Hi j0e*p(- (Kid  ♦ Xio)t) 
Ni3  (t)  * Nj30exp(-Kid3  t) 


The  differential  equations  form  the  mathematical  model  that 
was  tested  using  the  response  equations  in  Table  2.  The  viable 
spores  that  have  not  been  activated  (N3)  can  not  be  measured 
directly.  Activated  spores,  N2,  and  the  less  heat-r-esistant 
population,  N3  produce  colonies  in  an  enumeration  medium;  hence, 
they  are  the  survivors  of  the  heat  treatment.  Therefore,  the 
number  of  survivors,  N,  equals  the  summation  of  II2  ano  H3. 

.(31) 


N2  * N3 


In  consequence, 


equation  for 


survivors  (N) 


is  a summation  of  three  decaying  exponential  terms,  one  of  which 
is  negative: 

H(t)=(NlO*N2o)exp(-Kdt)-hioexpi-(KdrKo)t)*N3uexp(-Kd3t)..(J2) 

The  term  that  accounts  for  the  less  heat-resistant  fraction 
produces  an  initial  drop  in  the  number  of  survivors.  The  other 
positive  term  pertains  to  the  traditional  single  exponential 
describing  the  predominant  inactivation.  This  exponential  term 
takes  longer  to  decay  than  the  other  two,  and  it  becomes 
dominant  after  an  initial  period.  The  activation  part  of  tne 
survivor  curve  Is  character  iced  by  an  initial  rise  followed  by 
a transition  zone  (sometimes  referred  to  as  "shoulder”)  and 
later  by  a decrement  in  the  number  of  survivor  spores  (see 
Figure  4).  It  Is  determined  by  tne  combination  of  the  two  terms 
that  account  for  the  activated  spores.  The  negative  exponential 
explains  the  initial  rise  in  the  number  of  survivors.  There- 
fore, the  mathematical  model  should  be  able  to  describe  changes 

treatment,  and  the  physical  meaning  of  the  exponential  terms  in 


response  equations  is  clear. 


Experimental  Apparatus 


The  experimental  apparatus  permitted  experiments  for 
parameter  determination  and  validation  of  the  model.  It 
consisted  of  a three-neck  flask  equipped  with  a needle 
thermocouple,  a reflux  condenser  and  a set  of  9 needles  used  to 
inject  an  aliquot  of  the  spore  suspension  and  to  retrieve 
samples,  as  depicted  in  Figure  10.  The  thermocouple  was 
connected  to  an  analog-to-aigital  (a/d)  interface  By  means  of  an 
amplifier  board.  The  a/d  interface,  in  time,  was  connected  to 
an  IBH-pc  compatible  microcomputer  and  to  a heater  placed  under 
the  three-neck  flask.  This  computer  control  system  permittee 
accurate  control  of  the  temperature  of  the  dilute  spore 
suspension  inside  the  flask.  The  suspension  was  continuously 
mixed  with  a magnetic  stirrer,  while  time-temperature  data  were 
recorded  in  a computer  file.  A condenser  placed  in  tne  middle 
neck  recovered  evaporated  water  and  maintained  a constant  volume 
of  water  in  the  flask. 


Procedure 

The  flask,  the  magnet,  needles,  condenser,  tnermocouple,  99 
ml  of  distilled  water  and  a syringe  for  injecting  the  spore 
suspension  were  sterilized  before  each  experiment.  Sterile, 
disposable  syringes  were  used  for  sampling.  Special  care  was 


Figure  10.  Experimental  apparatus:  Three-ne 
flask  equipped  with  a needle 
thermocouple,  a reflux  condenser 
a set  of  nine  sampling  needles. 
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taken  to  assure  mixing  of  tne  spore  suspension  by  using  a vortex 
genie  for  at  least  20  s before  taxing  each  sample. 

The  apparatus  was  allowed  to  reach  equilibrium  at  the  desired 
treatment  temperature,  and  then  a 1.0  ml  aliquot  of  spore 
suspension  at  room  temperature  was  injected  into  the  99  ml  of 
sterile,  distilled  water  through  one  of  the  needles  installed  in 
the  left  neck  of  the  flask.  This  changed  the  internal 
temperature  of  the  flask  less  than  U,  and  the  process 
temperature  was  regained  after  a few  seconds.  No  attempt  was 
made  to  measure  or  regulate  the  pH  (pH  of  distilled  water  is  5 - 
5.5;  the  heat  resistance  of  spores  at  ph  7 is  expected  to  be 
higher).  Eignt  successive  0.5  ml  samples  were  taken  in  each 
experiment,  and  the  corresponding  times  were  recorded.  A 
different  needle  was  used  for  each  sample.  Similar  samples  were 
taken  of  tne  untreated  suspension  as  a control.  The  samples 
were  immediately  mixed  with  4.5  ml.  of  cold  sterile  water,  and 
kept  on  ice  or  in  a cold  room  until  the  microbiological 
enumeration  of  the  survivors  could  be  performed.  To  evaluate 
the  dispersion  of  the  data,  four  repetitions  of  one  point  were 
made  in  order  to  provide  an  estimate  of  the  standard  deviation 
used  later  in  the  nonlinear  regression  analysis. 

To  obtain  Initial,  tentative,  values  for  parameters  in  the 
model,  as  required  for  nonlinear  regression,  experimental  data 
were  analyzed  by  the  method  of  successive  residuals  (Hohsenin, 
1978)  which  is  very  useful  for  curve  fitting  when  the  model  is  a 


summation  of  exponentials. 


difference  in  the  exponential 


time  constants  determines  their  relative  rate  of  decay.  After  a 
sufficient  interval  the  slowest  exponential  in  the  sum 
dominates,  and  it  can  be  determined  by  regression.  When  this 
exponential  term  has  been  determined,  the  experimental  values 
can  be  subtracted  from  it,  leaving  another  exponential  to 
dominate.  It,  too.  can  be  determined  by  regression.  This 
procedure  is  applied  successively  until  the  differences  between 
the  experimental  values  and  the  last  exponential  determined 
become  negligible. 

The  procedure  to  find  good  inicial  values  for  tne  nonlinear 
regression  is  illustrated  below  for  a two-exponential  case 
(activation-inactivation).  The  equation  is 

y • aj  exp(-bj  t)  - a2  exp(-b;  t)  (33) 

where  y is  the  number  of  survivor  spores  at  time  t.  The  region 
where  the  positive  exponential  became  dominant  is  found  by 
inspection  of  the  survivor  curve,  and  the  values  of  aj  and  b[ 
are  determined  by  linear  regression  of  the  natural  logarithm  of 
N vs.  time  in  that  region.  Since  the  positive  exponential  and 
the  values  of  N are  known,  the  negative  exponential  can  be 
determined  by  subtraction: 


l-b2  t)  = ai  exp(-bi  t)  - 


.(34) 
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This  second  exponential  Is  referred  to  as  the  first  residual. 

The  method  of  successive  residuals  is  ver,  sensitive  to  changes 
in  slope.  Dispersion  of  the  experimental  data  In  this  study  was 
too  large  for  exclusive  use  of  tnis  method,  but  it  was  utilized 
to  find  good  first  estimates  for  use  In  the  nonlinear  regression 


procedure. 

A computer  program  was  developed  for  nonlinear  regression 
analysis  of  the  data  using  the  Levenberg-Harquardt  method, 
following  procedures  described  by  Press  et  al.,  (1986).  The 
value  of  chi  squared  was  used  as  the  criterion  of  the 
optimization  procedure  that  searched  for  a minimum,  and 
convergence  was  considered  to  be  achieved  when  the  change  in 
successive  iterations  was  smaller  than  10  ®.  In  order  to  reduce 
the  numerical  error  (round  off),  time  in  the  experiments  was 
divided  by  ten,  and  experimental  numbers  of  survivors  were 
divided  by  one  million.  This  scaling  of  variables  was  required 
to  facilitate  use  of  the  computer  software  and  to  improve  its 
accuracy  and  stability. 

Isothermal  experiments  were  performed  at  87,  93,  95,  and  99 
deg.  C Values  obtained  for  the  rate  constants  were  fitted  to  an 
Arrhenius  model  to  determine  the  dependence  of  rate  constants  on 
temperature.  The  Arrhenius  type  formulae  were  combined  with  the 
response  equations  to  simulate  nonisothermal  cases.  The 
nonisothermal  experiment  was  similar  to  the  isothermal 
experiment,  except  that  the  temperature  was  deliberately  varied 


with  time  by  turning  the  heater  on  and  off  several  times  during 
a process  while  keeping  an  accurate  record  of  the  time- 
temperature  history  thus  experienced. 

In  aouition  to  tne  data  generated  from  this  experimental 
work,  another  set  of  data  was  used  to  test  the  model.  These  data 
for  Clostridium  botulinum  in  oil  were  obtained  from  Dr.  L. 
Ababouch  at  the  Department  de  Microbiologie  Alimentalre  et 
Biotecnologle  of  the  Instltut  Agronomlque  et  Veterinaire  Hassan 
II  In  Morocco. 

The  theoretical  and  methodological  basis  of  this  research  led 
to  the  design  and  execution  of  a series  of  experiments.  The 
results  obtained  are  presented  and  discussed  in  the  next 


RESULTS  AND  DISCUSSION 
Experiments  with  Bloindicator  Units 

Experiments  performed  in  the  first  part  of  this  research  used 
a suspension  of  spores  of  the  thermophile  Bacillus  stearo- 
thermophilus  as  the  bloindicator  organism  inside  cylindrical 
bloindicator  units  (SlUs) . The  BIUs  were  made  of  plastic  or 
aluminum.  These  experiments  explored  the  importance  of  heat 
transfer  to  using  BIUs  in  the  validation  of  thermal  steriliza- 
tion processes.  Results  obtained  are  presented  in  the  following 

Heat  Transfer 


Time-temperature  graphs  shown  in  Figure  II  represent  typical 
temperatures  inside  the  BIUs  and  the  retort  used  in  these 
experiments.  There  was  a point  in  time  when  the  temperature 
inside  a plastic  BIU  rose  above  the  retort  temperature.  This 
happened  when  pressurized  air  was  allowed  to  enter  the  retort  at 
the  beginning  of  the  cooling  phase.  The  rate  of  heating  of 
plastic  containers  was  always  less  than  that  of  aluminum  con- 
tainers. The  ability  of  BIUs  to  follow  the  thermal  history  of 


Figure  11.  Time-temperature  curves  for  toe  retort; 
plastic  and  aluminum  Bioindicator 
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Figure  12. 


Heat  penetration  curves  for  the  plastic 
BIU  and  the  contents  of  the  can  using  an 
inverted  semilogarithmic  form  (Ball,  1957 
Stumbo,  1965) 


TEMPERATURE  (°F) 


TIME  (MINUTES) 


Figure  13. 


Heat  penetration  curves  for 
8IU  and  the  contents  of  the 
inverted  semilogarithmic  for 
St umbo,  1965) 


TEMPERATURE  (°F) 


TIME  (MINUTES) 


comparison  of  temperatures  Inside  a bioindicator  unit  and  the 
corresponding  point  in  a can  is  presented  in  the  form  of 
standard  heat  penetration  curves.  This  presentation  emphasizes 
tne  differences  that  may  occur  wnen  the  product  temperature  is 
near  the  retort  temperature.  Contributions  to  process  lethality 
in  this  region  are  very  important  (Stumbo,  1965,  page  122). 
Figure  13  shows  that  the  temperature  experienced  Inside  an 
aluminum  8IU  followed  the  external  temperature  history  very 
closely  (the  curves  are  nearly  identical);  whereas,  at  some 
points  in  time  within  the  lethal  range,  the  temperature  inside  a 
plastic  BID  appreciably  lagged  the  product's  temperature. 
Aluminum  appeared  to  have  no  detrimental  effect  on  the  spore 
suspension,  since  plate  counts  following  a given  thermal  process 
were  of  the  same  order  of  magnitude  in  both  plastic  and  aluminum 
Bills.  The  plate  count  information  is  presented  later. 

The  unexpected  rise  in  internal  temperature  of  a plastic  unit 
coincided  with  the  time  at  which  high  pressure  air  was  admitted 
to  the  retort  to  avoid  expansion  of  the  can  just  prior  to 
admitting  cooling  water.  In  order  to  confirm  this  observation, 
additional  experiments  were  conducted  with  tne  Bills  exposed 
directly  to  the  retort  environment  under  various  cooling  cycles 
witn  and  without  over-riding  air  pressure.  These  experiments 
always  showed  a rise  in  temperature  in  a plastic  bioindicator 
only  when  over-riding  air  pressure  was  admitted  to  the  retort. 


temperature  in  a plastic 


Figure  14.  Internal  temperatures  of  plastic  a 
aluminum  BIUs  when  no  pressurised 
is  admitted  into  the  retort. 


103 


rise  when  pressurized  air  was  not  used.  A rise  in  temperature 
was  not  observed  with  aluminum  containers  when  over-riding  air 
pressure  was  admitted.  These  results  prompted  an  examination  of 
the  mechanical  strength  properties  of  plastic  Bllis. 

Figure  15  presents  the  effect  of  temperature  on  the  force- 
deformation  curves  for  samples  prepared  from  plastic  BIOS.  The 
temperatures  used  were  lower  than  typical  temperatures  in  ster- 
. ilization  experiments  because  of  limitations  imposed  by  equip- 
ment. The  peak  value  decreased  from  1950  N at  70  deg.  C to  1270 
at  93  deg.  C.  When  a container  is  squeezed  by  external  pressure 
and  its  volume  is  reduced,  the  internal  temperature  should  rise 
as  is  explained  in  the  next  section  of  tnis  chapter. 

Process  Lethality 

The  values  for  0 at  250  deg.  F (6.73  min)  and  for  Z (18  deg. 
F)  were  determined  Independently  for  the  same  spore  suspension 
at  the  University  of  hassachusetts  (Evans,  1S84).  The 
lethalities  were  calculated  for  these  experiments  by  numerical 
integration  using  tne  variable  temperature  history  for  the  BIU 
investigation.  Lethalities  calculated  from  the  microbiological 
enumeration  of  survivors  and  from  integration  of  the  time- 
temperature  aata  are  summarized  in  Table  3: 


«•  £!s.sr.!ss?s;i«  iss 

from  plastic  Bills. 


105 


DEFORMATION  IN  *. 

• 70  OEC.  C.  ♦ 93  OEC.  C. 


Use 

Microbiological  Calculated 

KJSSSB” 

Kir""" 

8.0x10®  9.9x10® 

6.0x103  7.5x10® 

ssaisi..*-, 

2.5x10®  8.9x104 

“iSSSlEESlISKSafflS 


108 


understanding  the  problem.  Initially,  some  information  was 
available  from  the  literature  and  private  industry  on  the 
disagreements  between  final  spore  counts  found' using  plastic 
BIUs  containing  a B.  stearothermopnilus  spore  suspension  and  the 
number  of  surviving  spores  determined  by  simulation  with 
Bigelow's  model  for  the  thermal  history  of  a sterilization 
process.  A set  of  experiments  was  performed  to  verify  the 
accuracy  of  those  reports,  using  constant  and  variable 
temperature  processes.  The  discrepancies  in  survivor  spore 
counts  were  by  two  orders  of  magnitude  regardless  of  tne  type  of 
thermal  process  involved  ( see  upper  part  of  Table  4).  The 
performance  of  Bigelow's  model  was  not  good  enough  for  routine 
use  in  the  control  of  thermal  sterilization  processes. 

A second  set  of  experiments  at  constant  retort  temperature 
was  designed  to  test  the  influence  of  heat  transfer  on  the 
performance  of  BIUs.  Aluminum  BIUs  were  designed,  fabricated 
and  tested  against  plastic  BIUs.  The  experiments  yielded  time- 
temperature  histories  (figures  11,  12  and  13).  Figures  12  and 
13  are  presented  in  the  form  of  standard  industrial  heat 
penetration  curves  (Stumbo,  1965  ),  that  emphasize  tne 
temperature  region  where  the  internal  8IU  temperature  is  near 
the  retort  temperature.  This  region  is  where  a very  high  lethal 
effect  is  achieved.  Inspection  of  the  time-temperature  data 
confirmed  expectations  that  aluminum  BIUs  followed  tne  tnermal 
history  of  the  contents  of  a can  closely  enough  to  eliminate 
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heat-transfer  as  a source  of  discrepancy  in  spores  counts.  In 
this  set  of  experiments,  spores  were  not  heut-activated 
previous  to  their  exposure,  and  a vortex  genie  was  used  before 
and  after  exposure  to  stir  the  spore  suspension.  The  numbers  of 
surviving  spores  in  these  experiments  are  presented  in  the 
lower  part  of  Table  4.  The  discrepancies  in  spore  counts  are 
between  one  and  two  orders  of  magnitude,  and  they  are  far  larger 
than  the  15X  variation  among  replicates.  Inspection  of  thermal 
histories  revealed  an  interval  of  pressure  increases  at  the 
onset  of  cooling  that  caused  the  internal  temperature  of  the 
plastic  BIUs  to  rise  momentarily  above  the  constant  retort 
temperature  (250  deg.  F).  This  led  to  tne  design  and  execution 
of  experiments  where  this  phenomenon  was  investigated.  The 
temperature  rise  does  not  occur  when  pressure  is  not  Increased 
(Figure  14).  Temperature  induced  changes  in  the  mechanical 
strength  of  plastic  Bills  were  investigated  up  to  93  deg.  C.  The 
shape  of  force-deformation  curves  show  an  expected  softening  of 
the  material  at  high  temperatures.  The  contents  of  the  Bids 
included  some  air,  and  the  temperature  of  the  gaseous  phase  is 
expected  to  rise  when  the  volume  of  the  air-water  vapor  mixture 
is  reduced  due  to  changes  in  pressure.  To  evaluate  this  possi- 
bility, a gas  (air)  under  adiabatic  compression  and  a saturated 
air-water  mixture  were  analyzed. 
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Adiabatic  compression  is  described  by  tbe  following  equation 
(Vennard  and  Strett,  1982): 

Pl/UPl/RTi)*)  - P2/((P2/»T2)k)  (35) 

Pi  and  P2  are  the  absolute  internal  pressures  before 
and  after  squeezing,  respectively; 

Ti  and  T2  are  the  corresponding  absolute  temperatures; 

R is  the  universal  gas  constant;  and 

k is  the  adiabatic  exponent  for  air  (1.4). 

Calculations  with  this  equation  showed  that  a temperature 
rise  of  6 deg.  F (3  deg.  C)  could  be  expected  for  one  psi  in- 
crease in  internal  pressure.  Such  a pressure  increase  logically 
could  be  expected  from  the  facts  that  all  retort  controls  were 
operated  manually  and  the  mechanical  resistances  of  both  the  can 
and  the  Bill  would  reduce  volume  changes  to  some  extent.  A 
computer  program  (Sat-line.pas  in  Appendix  "A")  was  developed  to 
calculate  the  corresponding  changes  in  temperature  due  to  a 
pressure  rise  for  saturated  air  (air-water  vapor  mixture).  The 
results  in  this  case  showed  an  increment  of  3 deg.  F per  psi. 
(0.000242  deg.  C per  Pa)  (formula  from  ASAE,  1984). 

The  significance  of  these  findings  is  that  spore  suspensions 
in  plastic  Bills  under  the  pressure  conditions  previously 
described  could  experience  greater  lethality  than  in  the 


processed  product.  This  possibility  would  be  minimized  with 
rigid  aluminum  bioindicators.  Tnese  results  also  appear  to 
offer  a direct  corollary  to  the  work  reported  by  Teixeira  et 
al . , (1985).  and  Cleland  and  Gersterkamp  ( 1983)  who  investigated 
the  rapid  cooling  effect  of  sudden  retort  depressurization. 
Sudden  drops  in  retort  pressure  could  cause  less  tnan  tne 
expected  lethality  to  be  achieved  in  a process. 

The  experiments  performed  using  plastic  and  aluminum  BIUs 
were  very  significant  because  a)  the  relative  neat  transfer 
properties  could  not  explain  the  persistent  discrepancies  be- 
tween microbiological  and  calculated  values  of  the  number  of 
survivor  spores  and  bj  the  aluminum  8IUs  closely  followed  the 
thermal  history  of  the  contents  of  a can  and  provided  mechanical 
protection  in  case  of  a sudden  pressure  rise. 

Experiments  on  Kinetics 

Parameters  involved  in  the  response  formulae  (Table  2)  were 
determined,  and  tne  efficacy  of  the  response  formulae  in 
fitting  experimental  data  was  tested  using  a series  of 
isotnermal  experiments.  Subsequently,  the  model  was  validated 
by  comparison  of  predicted  and  experimental  responses  for  a 
nonisothermal  operation  of  the  sterilization  system. 


Isothermal  Experiments  for  parameter  determination 


The  number  of  surviving  activated  spores  as  a function  of 
time  for  treatments  at  four  constant  temperatures  are  presented 
In  Figures  16,  17,  IS  and  19.  Because  several  points  at  the 
beginning  of  tne  experiments  were  hidden  by  the  vertical  axis, 
the  origin  of  the  time  axis  has  been  slightly  displaced  in  these 
figures.  The  temperatures  of  the  treatments  were  87,  93,  95, 
and  99  deg.  C.  For  each  temperature,  two  curves  are  presented: 
one  represents  the  survivors  found  by  incubation  at  a tempera- 
ture of  45  deg.  C,  and  the  other  corresponds  to  survivors  found 
by  incubation  at  32  deg.  C.  The  general  trend  is  for  tne  two 
curves  to  begin  to  look  alike  at  high  temperatures  and  for 
spores  enumerated  at  32  deg.  C,  initially  constituting  a very 
small  number,  to  become  the  dominant  form  over  time.  With  the 
exception  of  the  treatment  at  87  deg.  C,  when  enough  time  has 
elapsed,  the  number  of  injured  survivor  spores  is  larger  than 


best  be  seen  by  viewing  all  four  cases  simultaneously,  as  in 
Figure  20. 

For  all  experiments,  activation  dominated  tne  initial  stages 
of  heat  treatment,  due  to  large  initial  increments  in  the  number 
of  survivors.  However,  at  higher  temperatures,  samples  could 
not  be  taken  sufficiently  fast  to  demonstrate  the  rapid 
increments  in  the  number  of  survivors  because  of  tne  logistics 


16.  Survival  curves  for  the  treatment 
of  spores  of  8.  subtilis  at 
87  deg.  C.  including  the  curve  for 
injured  spores  that  are  evident  by 
incubation  at  32  deg.  C and  normal 
spores  incubated  at  45  deg.  C. 
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Figure  17.  Survival  curves  fc 


Survival  curves  for  the  treatment 
of  spores  of  8.  subtills  at 
95  deg.  C,  including  the  curve  for 
injured  spores  that  are  evident  by 
incubation  at  32  deg.  C and  normal 
spores  incubated  at  45  deg.  C. 
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Figure  19.  Survival  curves  for  the  treatment 
of  spores  of  8.  subt i 1 is  at 
99  deg.  C,  including  the  curve  for 
injured  spores  that  are  evident  by 
incubation  at  32  deg.  C and  normal 
spores  incubated  at  <5  deg.  C. 
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(Thousands) 

Time  in  seconds. 

Incubation  at  32  C.  4 Incubation  at  45  C. 


Figure  20.  The  four 
together . 


c)  95  deg 


sets  of  survival  curves 
The  treatment  temperatures 
deg.  C,  b)  93  deg.  C, 
i.  C,  d)  99  deg.  C. 
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of  the  experimental  procedure.  It  follows  that  accuracy  of 
measurement  during  the  first  seconds  of  each  experiment  was 
limited.  Consequently,  the  possibility  of  an  initial  maximum 
occurring  near  the  beginning  of  a test  such  as  can  be  seen  in 
Figure  18,  could  not  be  confirmed  at  the  higher  temperature 
(figure  19).  In  summary,  the  initial  rise  is  determined  by 
activation,  and  it  becomes  snorter  at  higher  treatment  tempera- 


The  exponential  term  in  tne  response  formula  that  accounts 
for  the  less  heat-reslsteut  population  (Ng)  of  the  spore 
suspension  was  not  significant  in  this  research  because  the 
spore  suspension  used  was  homogenous  for  the  following  two 
reasons.  One,  the  spore  suspension  used  was  more  than  three 
years  old,  free  of  vegetative  cells  and  phase  daric  spores.  Two, 
the  organism  tested  did  not  present  multiple  genomes  (Lindsay, 

The  experimental  data  were  analyzed  using  nonlinear  regres- 
sion to  determine  values  of  the  parameters  in  equation  27. 
Calculated  values  for  the  survivor  spores  are  presented  together 
with  corresponding  experimental  data  in  figures  21,  22,  23, 

24,  25,  26,  27  and  28.  The  graphs  were  prepared  to  emphasize 
comparison  of  calculated  values  and  experimental  data.  Note 
that  the  scales  are  different  in  each  case.  The  numbers  pre- 
sented were  modified  to  minimize  round-off  error  during  the 

time  is  divided  by  10.  and  tne  number  of  survivors 


calculations: 


Comparison  between  calculated  survivor 
curves  and  experimental  data  for  spores 
of  B.  subti  l is  exposed  to  beat  treatment 
at  31  deg.  C and  incubated  at  45  deg.  C. 
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EXPERIMENTAL 


TIME 

CALCULATED 


Figure  22.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  injured 
8.  subtilis  spores  exposed  to  beat 
treatment  at  87  deg.  C. 
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Figure  23.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  spores 
• subtil  is  exposed  to  heat  treatment 


t §3  deg.  C and  incubated  ai 


EXPERIMENTAL 


TIME 

- CALCULATED 


Figure  24.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  injured 
B.  subtilis  spores  exposed  to  heat 
treatment  at  93  deg.  C. 


131 


TIME 

■ EXPERIMENTAL  CALCULATED 


Figure  25.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  spores 
of  8.  subtills  exposed  to  heat  treatment 
at  95  deg.  C and  incubated  at  45  deg.  C. 
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Figure  28.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  injured 
B.  subtilis  spores  exposed  to  heat 
treatment  at  95  deg.  C. 
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27.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  spores 
of  8,  subtilis  exposed  to  heat  treatment 
at  99  deg.  C and  incubated  at  45  deg.  C. 


TIME 

■ EXPERIMENTAL  CALCULATED 


Figure  28.  Comparison  between  calculated  survivor 
curves  and  experimental  data  for  injured 
6.  subtilis  spores  exposed  to  beat 
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by  one  million.  A natural  logarithmic  presentation  was  chosen 
in  agreement  with  the  exponential  nature  of  the  mathematical 


nonlinear  Curve-Fitting 

A computer  program  (Nonlin.pas  in  Appendix  "a")  was  developed 
for  nonlinear  least-square  curve  fitting  of  the  data.  Initial 
parameters  were  determined  from  samples  of  the  spore  suspension 
taken  at  the  beginning  of  each  experiment  and  arithmetic  means 
were  used  for  the  calculations:  5.9x10?  for  the  spores  incubated 
at  45  deg.  C and  1.3x10s  for  the  injured  spores;  the  standard 
deviations  were  1.72x10?,  and  1.77x10s  correspondingly. 

The  error  (standard  deviation  expressed  as  a percentage  of  the 
mean  value)  was  estimated  from  four  replicas  of  one  point  on  tne 
survivor  curve  to  be  201  and  30t  of  the  mean  for  the  injured 
spores.  Estimation  of  the  error  is  important  to  the  nonlinear 
regression  because  the  parameter  minimized  in  the  procedure  is 
Chi  squared,  and  a good  estimation  of  the  error  facilitates  the 
curve  fitting  process.  Increasing  the  standard  deviation  did 
not  affect  the  values  of  the  constants  determined  up  to  a loot 
of  the  value  of  the  mean.  Also,  the  rate  of  convergence 
subjectively  estimated  by  the  computation  time  did  not  change 
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significantly.  Low  susceptibility  of  the  nonlinear  curve 
fitting  process  to  error  allowed  use  of  a few  experimental 

The  same  procedure  applied  to  data  from  Clostridium  Qotulinum 
62A  in  olive  oil  provided  by  Or.  L.  Ababouch  tested  the  part  of 
the  model  that  accounts  for  a less  heat-resistant  population. 

The  results  are  presented  in  Figure  29,  and  the  good 
correspondence  between  the  values  calculated  with  the  response 
formulae  and  the  experimental  data  is  obvious.  In  this  case, 
the  treatment  temperature  was  120  deg.  C 

A scatter  diagram  of  all  corresponding  calculated  and 
experimental  values  is  presented  in  Figure  30.  The  straight 
line  represents  the  theoretical  case  of  perfect  correlation 
between  predicted  and  actual  responses.  Tne  regression  equation 
(not  presented  in  the  graph)  is 

In  y = 0.964683  In  x * 0.0303  (36) 

where  y and  x represent  calculated  and  experimental  values,  and 
the  correlation  coefficient  (r)  is  0.96343.  with  76  degrees  of 
freedom.  The  slope  very  close  to  1 and  the  very  small  intercept 
confirm  the  efficacy  of  the  proposed  moael  in  describing  the 
dynamics  of  spore  populations  under  constant  lethal  temperature 
treatment.  A region  of  some  discrepancy  is  found  in  the  initial 
moments,  as  can  be  seen  in  figures  21  and  25;  this  region  should 


Comparison  between  calculated  values 
and  experimental  data  for 
Clostridium  botulinum  62  A in  olive  oil. 
(Experimental  data  provided  by 
L.  Ababouch,  1986). 
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EXPERIMENTAL 


CALCULATED 


Scatter  diagram  showing 
correlation  between  experimental 
and  calculated  values  for  a total 
of  68  separate  sets  of  experimental 
conditions. 
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(perimQntal  values. 
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be  the  subject  of  additional  research.  High  temperature-short 
time  processes  may  benefit  from  a better  understanding  of  the 
events  that  occur  during  the  first  crucial  moments  of  exposure 
to  lethal  temperature. 

Temperature  dependence  of  rate  constants.  The  rate  constants 
of  activation  and  Inactivation  for  the  standard  case  and  for 
spores  incubated  at  32  deg.  C are  presented  in  figures  31,  32, 
33,  and  34.  The  graphs  follow  the  form  of  the  Arrhenius' 
equation  for  temperature  dependence  of  rate  constants,  where  the 
natural  logarithm  of  the  rate  constant  is  plotted  vs.  the  recip- 
rocal of  the  absolute  temperature  (1/T).  Activation  energy  (Ea) 
was  calculated  from  the  slope  of  the  regression  equation  for  a 
graph  of  the  natural  logarithm  of  a rate  constant  as  a function 
of  the  reciprocal  of  absolute  temperature.  The  regression 
equations  are  as  follows: 

For  the  spores  incubated  at  45  deg.  C, 

In  Kq  » - 27201. 2/T  ♦ 65.64113;  r = - 0.9989  (37) 

In  Ka  ■ - 17808. 9/T  ♦ 41.50050;  r = - 0.9278  (38) 

For  the  spores  incubated  at  32  deg.  C, 

In  Kj<j  = - 34529.4/T  ♦ 84.90411;  r = - 0.9974  ...(39) 

In  Kj  ■-  43547. 4/T  ♦ 113.1216;  r •-  0.9352  ...(40) 

The  degrees  of  freedom  are  two  in  all  cases. 


Arrhenius  plot  for  the  inactivation 
rate  constant  for  spores  of  8.  subtilis 
incubated  at  45  deg.  C. 
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1/ (ABSOLUTE  TEMPERATURE)  . 
REGRESSION  LINE. 


Figure  32.  Arrhenius  plot  for  the  inactivation 
rate  constant  of  B.  subtilis  injured 
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REGRESSION  LINE. 


Figure  33.  Arrhenius  plot  fi 


Activation 
ss  of  B-  subtllis 
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1 / (ABSOLUTE  TEMPERATURE) . 
REGRESSION  LINE. 


Figure  34.  Arrhenius  plot  for  the  activation 

rate  constant  of  B.  subtiiis  injured 
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1 / (ABSOLUTE  TEMPERATURE) . 


REGRESSION  LINE. 


The  Arrhenius  equation  may  be  written  as  follows  ( Aiba  and 


Humphrey,  1973) 

K ■ a'  expi-  Ea/(RT));  or 

In  K = In  a'  - Ea/(RT)  (41) 

The  slope  (m)  of  the  curve  is  then 
m = - Ea/«  (42) 


The  value  of  the  universal  gas  constant  ,R,  is  8314 
J/(kmol  K)  (Loncin,  1977). 

The  activation  energy  is  then 


Ea  = - mR  .(43) 

Values  determined  for  the  activation  energy,  in  J/kmol,  are 
presented  in  Table  5. 


Table  5.  Activation  energy  (J/kmol). 


Inactivation 

Activation* 

Inactivation** 

Activation** 


"incubation  at  32  deg.  C 


2.26x11)8 

1.48x10* 

2.87x10* 

3.62x10* 


Dependence  of  the  rate  constant  of  inactivation  for  the  two 
incubation  temperatures  of  45  and  32  deg.  C is  described  very 
well  by  the  Arrhenius  equation,  and  the  dispersion  of  the  data 
is  very  low.  Activation  energy  corresponds  to  i in  the  Qxo 
model.  Determination  of  the  activation  constant  is  influenced  by 
logistical  limitations  of  the  experimental  procedure.  The  first 
two  points  in  figures  32  and  33,  corresponding  to  higher  treat- 
ment temperatures,  show  some  disagreement.  This  discrepancy 
emphasizes  the  need  to  explore  the  high  temperature-short  time 
region  in  order  to  provide  better  values  for  tne  constants  in 
the  Arrhenius  equation  describing  the  temperature-dependence  of 
activation.  Prediction  of  activation  rate  constants  was  not  as 
accurate  as  prediction  of  the  inactivation  rate  constants,  but 
the  absolute  value  of  the  correlation  coefficient  was  higher 


1 Validation— Nonisothermal  Experiment 


The  model  was  tested  by  ar 
io  longer  was  held  constant, 
ient  temperature  used  was  re 


experiment  wherein  the  temperature 
rhe  tnermal  history  of  the  tran- 
:orded  by  conputer  and  it  is 


presented  i 


The  results  support  usi 


traditional  model. 


Figure  35.  Time-temperature  history 
nonisothermal  experiment. 
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TIME  IN  SECONDS. 


gradually,  probably  due  to  the  removal  and  use  of  the  lower 
spore-rich  populations  In  the  tube  that  held  the  spore  suspen- 

microseopic  count  was  9.7  x 107  and  tne  standard  deviation  was 
2.3  x 10®. 

A computer  program  (Noniso.pas  in  Appendix  "C")  was  developed 
to  calculate  the  number  of  survivor  spores  that  result  from  a 
given  time-temperature  history  (simulation)  using  the  improved 
model.  The  rate  constants  were  calculated  from  the  regression 
equations  for  the  Arrhenius  formula,  and  the  initial  populations 
used  were  the  averages  shown  in  Table  6.  The  program  calculated 
changes  in  the  spore  population  for  each  time  increment,  and  it 
used  the  new  populations  obtained  as  initial  values  for  the  next 
time  interval  and  recorded  them  in  a computer  file,  together 
with  the  time  in  seconds. 

A comparison  of  the  number  of  surviving  spores  calculated 
from  the  response  formulae  and  the  values  determined  experiment- 
ally is  presented  in  Figure  36.  The  numbers  are  always  of  the 
same  order  of  magnitude,  and  there  are  positive  and  negative 
deviations  from  tne  ideal  case.  This  is  a significant  improve- 
ment over  the  results  obtained  using  the  traditional  Sigelow 

The  logarithmic  nature  of  lethality  makes  it  very  sensitive 
to  changes  in  the  value  of  Djjjg.  Values  of  D and  z used  in 
this  research  were  inadequate  to  predict  the  dynamics  of  a 


nonisothemal  experiment. 
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EXPERIMENTAL  DATA. 
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population  of  bacterial  spores  using  Bigelow's  model.  Improved 
kinetic  models  describing  inactivation  of  bacterial  spores  could 
play  a major  role  in  achieving  more  precise  prediction  of  spore 
survival  in  a lethal  heat  treatment. 

Discussion.  Development  of  the  mathematical  model  in  this 
study  was  supported  by  a series  of  isothermal  experiments,  using 
B.  subtilis  spores,  to  determine  values  of  the  parameters  Inv- 
olved and  to  verify  its  ability  to  respond  properly  to  constant 
temperature  treatment.  The  temperature  dependence  of  rate 
constants  also  was  determined,  and  a variable-temperature  exper- 
iment tested  the  capability  of  the  mathematical  model  with  the 
temperature  dependence  of  the  rate  constants  incorporated. 

Figure  20  presents  four  survivor  curves  together  in  order  to 
facilitate  comparison.  Logistical  limitations  of  tne  meas- 
urement technique  at  the  beginning  of  each  experiment  led  to 
some  uncertainty  in  the  activation  rate  constants,  particularly 
at  higher  temperatures  (Figures  32  and  33).  The  survivor  curve 
for  the  99  deg.  C treatment  may  have  reached  a maximum  value 
sooner  than  it  appears  in  the  graph.  This  maximum  would  be 
expected  as  a consequence  of  rate  constants  increasing  with 
temperature;  this  idea  is  supported  by  the  first  two  points  of 

Injured  spores  increase  in  relative  importance  with  time.  In 
general,  at  the  end  of  each  experiment  they  became  the  dominant 


population  of  survivors.  The  dominance  of  injury  in  the  long 
term  emphasizes  the  importance  of  including  it  in  a kinetic 
model.  Injury  may  be  particularly  important  in  cases  where 
storage  conditions  mignt  favor  the  germination  and  outgrowth  of 
injured  spores. 

The  capability  of  tne  model  to  describe  changes  in  the 
number  of  survivor  spores,  as  shown  in  Figure  30,  for  Cl.  botu- 
linum  spores  that  Included  a less-heat  resistant  population, 
emphasizes  the  potential  of  the  approach  followed  in  this  work 
to  handle  complex  situations  that  often  occur  in  food  engineer- 
ing and  the  importance  of  dealing  with  these  problems  in  tneir 
complexity.  The  final  test,  using  a variable  temperature  pro- 
cess, can  be  considered  successful  because  the  discrepancies 
found  can  be  attributed  to  experimental  error.  This  model 
should  be  used  in  order  to  achieve  meaningful  results  with 
microbiological  validation  of  thermal  sterilization  processes. 
This  improved  model  may  also  become  a tool  in  basic  research 
leading  to  an  even  more  realistic  mathematical  characterization 
of  bacterial  spore  behavior  under  treatment  at  lethal  tempera- 

Tne  need  for  an  improved  kinetic  mooel  is  supported  by  the 
existence  of  concomitant  reactions  as  reported  in  the 
literature.  Reactions  such  as  activation  and  injury  and  tne 
possivility  of  alternative  germination  mechanisms  and  less-heat 
resistant  populations  are  not  considered  in  Bigelow's  model. 
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Yet,  all  of  these  factors  definitely  Influence  the  spore  popula- 
tion surviving  a thermal  process.  Special  care  has  been  taken 
to  develop  computer  programs  that  permit  application  of  the  new 
model  by  microcomputer,  as  a means  of  facilitating  use  of  tne 


model  in  the  design  and  validation  of  thermal  sterilization 


CONCLUSIONS  ANO  SUGGESTIONS  FOR  FUTURE  WORK 
Conclusions 

The  aluminum  bioindicator  units  designed,  fabricated,  and 
evaluated  in  this  research  were  determined  to  be  superior  to 
plastic  units  used  in  commercial  practice  for  validating  thermal 
sterilization  process.  The  superiority  lies  in  the  superior 
heat  transfer  and  rigidity  of  the  aluminum  units.  However,  heat 
transfer  and  mechanics  can  not  explain  completely  the  large 
discrepancy  often  found  between  the  number  of  spores  surviving  a 
given  thermal  history  as  determined  by  Bigelow's  kinetic  model 
and  by  microbiological  enumeration. 

An  improved  mathematical  model  of  a bacterial  spore 
population  under  lethal  thermal  treatment  was  developed  and 
tested  satisfactorily.  The  new  model  took  Into  consideration 
three  concomitant  spore  transformations  (inactivation, 
activation  and  injury)  as  well  as  the  presence  of  less  heat- 
resistant  fractions  and  the  activation  of  alternative 
germination  mechanisms. 

The  understanding  of  processes  related  to  thermal 
inactivation  of  bacterial  spores  has  been  increased.  The 
following  statements  summarize  those  insights: 
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1.  Bigelow's  kinetic  model  is  significantly  imprecise  for 
validation  of  commercial  sterilization  processes  because  it  does 
not  reflect  the  relevant  biology. 

2.  Aluminum  bioindicator  units  follow  the  thermal  history  of 
the  contents  of  the  can  more  closely  than  plastic  units,  and  are 
less  susceptible  to  pressure-induced  temperature  changes. 

3.  System  analysis  approach  was  successfully  used  to  develop 
a model  of  a bacterial  spore  population  valid  in  the  lethal 
temperature  range. 

4.  The  transformations  studied  — activation.  Inactivation 
and  injury—  followed  pseudo-first  order  reaction  kinetics 
within  the  time  Intervals  used  in  the  experiments,  and  they 
behaved  as  Independent,  simultaneous  and  irreversible  reactions. 

5.  The  mathematical  model  developed  in  this  work  closely 
described  the  dynamics  of  the  number  of  survivor  spores  under 
treatment  at  constant  lethal  temperature. 

6.  Arrhenius'  equation  conveniently  describes  the 
temperature  dependency  of  the  rate  constants  in  the  temperature 

7.  Using  the  Arrhenius  equation  to  describe  the  temperature 
dependence  of  parameters  in  the  mathematical  model,  the  response 
formulae  predicted  changes  in  the  number  of  survivor  and  injured 
spores  when  the  temperature  was  allowed  to  vary  with  time. 

8.  For  the  spore  suspension  tested,  the  number  of  spores 
capable  of  growing  at  32  deg.  C was  higher  than  the  number  of 
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survivor  spores  determined  by  Incubation  at  45  deg.  C in  the 
final  range  of  the  time  interval  tested. 

9.  The  experimental  procedure  adopted  in  tnls  study  requires 
small  sample  sites,  is  simple  to  use,  and  yields  good  results 
for  temperatures  up  to  the  boiling  point  of  water.  It  also 
circumvents  many  of  the  shortcomings  of  traditional  TOT  cans  and 
glass  tubes  in  oil  or  steam  baths. 

10.  The  computer  programs  developeo  in  this  research  for 
data  analysis  and  simulation  show  the  great  potential  of  system 
analysis  and  micro-computers  as  tools  in  applied  microbiology. 

Suggestions  for  Future  Work 

The  order  of  reaction  should  be  studied  carefully  for  the 
transformations  considered  in  tne  model  presented  in  this 
dissertation.  Inspection  of  some  data  for  the  long  term  changes 
in  Bacillus  stearotnermophilus  (Lindsay,  1986)  shows  a departure 
from  first  order  kinetics  that  can  not  be  explained  by  a mixed 
population.  A more  accurate  description  of  the  kinetics  may 
provide  a better  mathematical  model.  The  following  are  specific 
suggestions: 

1.  The  early  moments  of  neat  treatment  at  lethal  temperature 
should  be  studied  in  more  detail.  8etter  experimental 
equipment,  allowing  faster  sampling  and  higher  treatment 
temperature,  can  be  developed  following  basic  concepts  of  the 
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apparatus  used  in  this  work:  computer  control  and  data 
acquisition,  injection  of  a small  sample  of  concentrated  spore 
suspension  into  an  already  equilibrated  and  well  mixed  reactor, 
and  handling  of  the  samples  at  low  temperatures.  The  main 
challenge  is  a very  fast  sampling  procedure  that  prevents 
contamination  across  large  pressure  differences. 

2.  The  procedure  followed  in  this  work  should  be  applied  to 
study  recovery  of  spores  from  injury  and  inactivation  related  to 
long-term  storage  of  canned  products. 

3.  In  general,  population  dynamics  and  system  analysis 
should  be  used  for  the  solution  of  complex  food  engineering 
problems. 
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APPENDIX  A 

COMPUTER  PROGRAM  TO  CALCULATE  SATURATION  TEMPERATURE 
SATURATION  PRESSURE  FOR  AIR-WATER  MIXTURES. 


:i,3TSr,sir,sxi"i,'r 
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(•file  SATLINE  .OAT  with  the  values  of  the  coefficients*) 


35.1579 

24.5926 

2.11821 

-0.341447 

-6.0313296 

0.00386583 

-0.249018E-3 

0.684016E-5 


APPENDIX  8 

COMPUTER  PROGRAM  FOR  THE  NONLINEAR  REGRESSION  ANALYSIS  OF 


DATA 


r 4-ft  r ft  || 


fiSss; to.. 


b( i row , 1 ) :=  b(icol,l); 
b(icol.l)  :=  dura 
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ind*r(i)  Irow; 
indxc(i)  :■  icol; 

IF  (a(icol.icol)  = 0.0)  THEN  BEGIN 

writeln( 'pause  2 in  GAUSSJ  - singular  matrix ' ) ; 
In 

ENO; 

piviny  :=  1.0/a(icol ,icol); 
a(icol.icol)  1.0; 

FOR  1 1 TO  n DO  BEGIN 

a(icol,l)  :=  a(icol,l)*pivinv 


FOR  11  :=  1 TO  n DO  BEGIN 
IF  (11  <>  icol)  THEN  BEGIN 
dum  :=  a(ll,icol); 
a(  1 1 .icol ) :■  0.0; 

FOR  1 :=  1 TO  n DO  BEGIN 

a(ll,l)  s»  a(ll,l)-a(icol,l)*dura 


k 1 TO  n 
ura  :=  a(k,indxr(l)); 
(k.indxrflj)  :=  a(k,indxc(l)); 


PROCEDURE  covsrt(VAR  covar;  glcovar;  ncvn:  integer;  ma: 
integer; 

lista:  gllista;  mfit:  integer); 

(•^Programs  using  routine  COVSRT  must  define  the  types 

glcovar  = ARRAY  (l..ncvm,l..nc*m)  of  real; 
gllista  = ARRAY  (l..mfit)  OF  integer; 
in  the  calling  program.  *) 


SiffiaK... 


j pss±"> 


"AwiiSJSS,  D?=*glatry(lista(j))*da(j) 


isar"!&«n...0,i 

""kTssr- 


ESK-S;;,,,™.,,  .alpha :glncabynca ;chisc 

”ir-  isijssfiin... 


J i:=l  TO  rnnadat 

DO  witeln(fout.,a(,,1,,J  = '.■{<):U:<}; 


ZOO 


> nriteln(fout); 

> writelnjfout, 'The  covariance  matrix  is:  ); 

> FOR  1 :=1  TO  mmadat  00 

> BEGIN 

> FOR  j:=l  TO  rnnadat 

> DO  write(fout,covar(i,j):l5:5,  ); 

> writeln(fout); 

> END;(*for  1*) 

> writeln(fout);  . 

> writelnjfout, 'The  curvature  matrix  is:  ); 

> FOR  i:=l  TO  rnnadat  DO 

> BEGIN 

> FOR  j:BI  TO  mmadat 

> 00  write(fout,alpha(i,j):15:4,'  ); 

> writeln(fout); 

> ENO;(*for  1*) 

> Hrlteln(fout);  . 

> writeln(fout,‘That  is  all  there  is,  there  am  t no  more 

>^’  close(fout); 

> ENO;(*reporte*j 


reset(fin); 

readln(f in, numdat .factdum) ; 
(•initialize  arrays*) 

FOR  i:-l  TO  numdat  DO 
BEGIN 

x(i):=0.0; 

y(1):-0.0; 

sig(i):=0.0; 

FOR  j:-l  TO  numdat  DO 
BEGIN 

covar(i,j):"0.0; 
alpha(i,j):=0.0; 
ENO;(*for  j*) 
EN0;(*for  i*) 

(•initialize  variables*)  ^ 


a(2):=0. 00122; 

a(3):=100.0; 

a(4):=0.1; 

mfitdat:=2; 

lista(l):=3; 

11sta(2):=4; 

(*  for  1:*1  tow 


procedures*) 

(•tentative  values  for  the 
parameters*) 


FOR  1:*  1 10  numdat 

M readln(fin,*(i).y(')»s’0(’))! 


(*  writeln(i,x(1));*) 

-®SKaM---s , 

" lista.mf'tdat ,covar,alpna,ni 


oldchisq:=chisq; 

criterio:=false; 

compare(oldchisq,chisq,cnterio). 


comparaiuiui-ii 

URITELN(CHISQ) ; 
oldchisq:=chisq; 
UNTIL  criterio: 

alaabda:=0.0; 

®rqmin(x. 


(•final  call  of 

mix  v.siu."anKJat,a,wiadat, 

lista  ,«f  i tdat  .covar  .alpha  ,r 

reporte(covar.alpha,chisq,a). 

close(fin); 

END.(*nonlin*) 


the  procedures*) 

:a,chisq,alambda); 
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